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Abstract 

Higgs boson production at the LHC within the Standard Model and its min- 
imal supersymmetric extension is reviewed. The predictions for decay rates 
and production cross sections are updated by choosing the present value of 
the top quark mass and recent parton density sets. Moreover, all relevant 
higher order corrections, some of which have been obtained only recently, are 
included in a consistent way. 
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1 Introduction 



1.1 Standard Model 

The Higgs mechanism is a cornerstone of the Standard Model (SM). To formulate the 
standard electroweak theory consistently, the introduction of the fundamental Higgs field 
is necessary |1|]. It allows the particles of the Standard Model to be weakly interacting up 
to high energies without violating the unitarity bounds of scattering amplitudes. The uni- 
tarity requirement determines the couplings of the Higgs particle to all the other particles. 
These basic ideas can be cast into an elegant and physically deep theory by formulating 
the electroweak theory as a spontaneously broken gauge theory. Due to the fact that 
the gauge symmetry, though hidden, is still preserved, the theory is renormalizable 
The massive gauge bosons and the fermions acquire their masses through the interaction 
with the Higgs field The minimal model requires the introduction of one weak isospin 
doublet leading, after the spontaneous symmetry breaking, to the existence of one ele- 
mentary scalar Higgs boson. Since all the couplings are predetermined, the properties of 
this particle are fixed by its mass, which is the only unknown parameter of the Standard 
Model Higgs sector. Once the Higgs mass will be known, all decay widths and production 
processes of the Higgs particle will be uniquely determined . The discovery of the Higgs 
particle will be the experimentum crucis for the standard formulation of the electroweak 
theory. 

Although the Higgs mass cannot be predicted in the Standard Model, there are several 
constraints that can be deduced from consistency conditions on the model Upper 
bounds can be derived from the requirement that the Standard Model can be extended up 
to a scale A, before perturbation theory breaks down and new non-perturbative phenom- 
ena dominate the predictions of the theory. If the SM is required to be weakly interacting 
up to the scale of grand unified theories (GUTs), which is of (9(10^^ GeV), the Higgs 
mass has to be less than ~ 200 GeV. For a minimal cut-off A ~ 1 TeV and the condition 
Mh < a, a universal upper bound of ~ 700 GeV can be obtained from renormalization 
group analyses 0, Q and lattice simulations of the SM Higgs sector 0. 

If the top quark mass is large, the Higgs potential may become unbounded from below, 
rendering the SM vacuum unstable and thus inconsistent. The negative contribution of 
the top quark, however, can be compensated by a positive contribution due to the Higgs 
self-interaction, which is proportional to the Higgs mass. Thus for a given top mass 
Mt = 175 GeV |^ H a lower bound of ~ 55 GeV can be obtained for the Higgs mass, 
if the SM remains weakly interacting up to scales A ~ 1 TeV. For A ~ Mqut this 
lower bound is enhanced to Mh ^ 130 GeV. However, the assumption that the vacuum 
is metastable, with a lifetime larger than the age of the Universe, decreases these lower 
bounds significantly for A ~ 1 TeV, but only slightly for A ~ Mqut 0- 

The direct search in the LEP experiments via the process e~^e~ —>■ Z*H yields a 
lower bound of ~ 77 GeV on the Higgs mass 0. This search is being extended at the 
present LEP2 experiments, which probe Higgs masses up to about 95 GeV via the Higgs- 
strahlung process Z* ZH p0|-p!2|. After LEP2 the search for the SM Higgs particle 
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will be continued at the LHC for Higgs masses up to the theoretical upper limit [|T3], . 
The dominant Higgs production mechanism at the LHC will be the gluon-fusion process 



15 



PP 99 



H. 



which provides the largest production cross section for the whole Higgs mass range of 
interest. For large Higgs masses the W and Z boson-fusion processes ]TB| , |T7| 



PP 



qq^qq + WW/ZZ Hqq 



become competitive. In the intermediate mass range Mz < Mh < 2Mz Higgs-strahlung 
off top quarks [1^ and W, Z gauge bosons [|T9|, |2y] provide alternative signatures for the 
Higgs boson search. 

The detection of the Higgs boson at the LHC will be divided into two mass regions: 



(i) 



For Mw < Mh ^ 140 GeV the only promising decay mode is the rare photonic one, 
H '-f'-f, which will be discriminated against the large QCD continuum background 
by means of excellent energy and angular resolutions of the detectors |14|. Alter- 



natively excellent /x-vertex detectors might allow the detection of the dominant bb 
decay mode pT[, although the overwhelming QCD background remains very diffi- 



cult to reject [0]. In order to reduce the background it may be helpful to tag the 
additional W boson in the Higgs-strahlung process pp —>■ HW |T9| , pO| or the tt pair 
in Higgs bremsstrahlung off top quarks, pp Hti 



(ii) 



In the mass range 140 GeV < < 800 GeV the search for the Higgs particle can be 
performed by looking for final states containing 4 charged leptons, which originate 



from the Higgs decay H 



Tl 



The QCD background will be small so 
that the signal can be extracted quite easily. For the Higgs mass region 155 GeV 
< Mh ^ 180 GeV another possibility arises from the Higgs decay H WW^*^ — > 
V 



23], because the W boson decay mode is dominating by more than one order 



of magnitude in this mass range, while the Z pair decay mode may be difficult to 
detect due to a strong dip in the branching ratio BR(if — > ZZ*) for Higgs masses 
around the W pair threshold. For Higgs masses above ~ 800 GeV the search may 
be extended by looking for the decay chains H — > ZZ, WW llvv. A Higgs boson 
search up to ~ 1 TeV seems to be feasible at the LHC ||14|. 



In order to investigate the Higgs search potential of the LHC, it is of vital importance 
to have rehable predictions for the production cross sections and decay widths of the Higgs 
boson. In the past higher order corrections have been evaluated for the most important 
processes. They are in general dominated by QCD corrections. The present level leads 
to a significantly improved and reliable determination of the signal processes involved in 
the Higgs boson search at the LHC. 
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1.2 Supersymmetric Extension 

Supersymmetric extensions of the SM p^, ^ are strongly motivated by the idea of pro- 



viding a solution of the hierarchy problem in the SM Higgs sector. They allow for a light 
Higgs particle in the context of GUTs ||26[, in contrast with the SM, where the extrapo- 



lation requires an unsatisfactory fine-tuning of the SM parameters. Supersymmetry is a 
symmetry between fermionic and bosonic degrees of freedom and thus the most general 
symmetry of the S-matrix. The minimal supersymmetric extension of the SM (MSSM) 
yields a prediction of the Weinberg angle in agreement with present experimental mea- 



surements in the context of GUTs ^7j. Moreover, it does not exhibit any quadratic 
divergences, in contrast with the SM Higgs sector. Throughout this review we will con- 
centrate on the MSSM only, although most of the results will also be qualitatively valid 
for non- minimal supersymmetric extensions ||28|| . 

In the MSSM two isospin Higgs doublets have to be introduced in order to preserve 
supersymmetry |2^. After the electroweak symmetry-breaking mechanism, three of the 



eight degrees of freedom are absorbed by the Z and W gauge bosons, leading to the 
existence of five elementary Higgs particles. These consist of two CP-even neutral (scalar) 
particles h, H, one CP-odd neutral (pseudoscalar) particle A, and two charged particles 
H^. In order to describe the MSSM Higgs sector one has to introduce four masses M^, 
Mh, Ma and Mh± and two additional parameters, which define the properties of the 
scalar particles and their interactions with gauge bosons and fermions: the mixing angle 
(3, related to the ratio of the two vacuum expectation values, tg/3 = V2/V1, and the mixing 
angle a in the neutral CP-even sector. Due to supersymmetry there are several relations 
among these parameters, and only two of them are independent. These relations lead to a 
hierarchical structure of the Higgs mass spectrum [in lowest order: Mh < Mz, Ma < Mh 
and Mw < Mh±]- This is, however, broken by radiative corrections, which are dominated 
by top-quark-induced contributions ^ . The parameter tg/5 will in general be assumed 
to be in the range 1 < tg/5 < mt/nib [vr/4 < /5 < n/2], consistent with the assumption 
that the MSSM is the low-energy limit of a supergravity model. 

The input parameters of the MSSM Higgs sector are generally chosen to be the mass 
Ma of the pseudoscalar Higgs boson and tg/5. All other masses and the mixing angle a 
can be derived from these basic parameters [and the top and squark masses, which enter 
through radiative corrections]. In the following qualitative discussion of the radiative 
corrections we shall neglect, for the sake of simplicity, non-leading effects due to non- zero 
values of the supersymmetric Higgs mass parameter /i and of the mixing parameters At 
and Ab in the soft symmetry-breaking interaction. The radiative corrections are then 
determined by the parameter e, which grows with the fourth power of the top quark mass 
Mf and logarithmically with the squark mass M5, 



v/27r2sin2/5 

These corrections are positive and they increase the mass of the light neutral Higgs boson 
h. The dependence of the upper limit of Mh on the top quark mass Mt can be expressed 
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as 

Ml < M| cos^ 2p + e sin^ p . (2) 

In this approximation, the upper bound on is shifted from the tree level value Mz up 
to ~ 140 GeV for Mt = 175 GeV. Taking Ma and tg/5 as the basic input parameters, the 
mass of the lightest scalar state h is given by 



Ml 



Mi +Mi 



{Ml + M| + e)2 - 4MlM| cos2 2/5 - 4e(M2 sin^ /5 + M| cos2 /?) 



(3) 



The masses of the heavy neutral and charged Higgs bosons are determined by the sum 
rules 



Ml + Ml-Ml + e 
Ml + M^. 



The mixing parameter a is fixed by tg/S and the Higgs mass Ma 

Ml + Ml 



tg2a = tg2/? 



Ml-M| + e/ cos 2/? 



with 



TT 



< a < . 



(4) 



(5) 



The couplings of the various neutral Higgs bosons to fermions and gauge bosons depend on 
the angles a and (3. Normalized to the SM Higgs couplings, they are listed in Table |l[ The 
pseudoscalar particle A does not couple to gauge bosons at tree level, and its couplings 
to down (up)-type fermions are (inversely) proportional to ig(3. 





«T> d) ii> 

9u 9d 9v 


SM H 


1 1 1 


MSSM h 
H 
A 


cos a/ sin /3 — sin a/ cos /3 sin(/5 — a) 
sin a/ sin /3 cos a/ cos /5 cos(/3 — a) 
1/tg/? tg(3 



Table 1: Higgs couplings in the MSSM to fermions and gauge bosons [V = W, Z] relative 
to SM couplings. 



Recently the radiative corrections to the MSSM Higgs sector have been calculated up 
to the two- loop level in the effective potential approach . The two-loop corrections are 



dominated by the QCD corrections to the top-quark-induced contributions. They decrease 
the upper bound on the light scalar Higgs mass M^ by about 10 GeV. The variation of 
Mh with the top quark mass is shown in Fig. ^ for Ms = 1 TeV and two representative 
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Figure 1: (a) The upper limit on the light scalar Higgs pole mass in the MSSM as a 
function of the top quark mass for two values of tgjS = 1.5,30. The top quark mass has 
been chosen as Mt = 175 GeV and the common squark mass as Ms = 1 TeV. The full 
lines correspond to the maximal mixing case [At — y/6Ms, — ji — and the dashed 
lines to vanishing mixing. The pole masses of the other Higgs bosons, H, A, H^, are 
shown as a function of the pseudoscalar mass in (b-d) for two values of tg(3 — 1.5, 30 and 
vanishing mixing. 
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values of tg/3 = 1.5 and 30. While the dashed curves correspond to the case of vanishing 
mixing parameters fi = At = Af, = 0, the solid lines correspond to the maximal mixing 
case, defined by the Higgs mass parameter /i = and the Yukawa parameters Ajj = 0, 
At = y/QMs. The upper bound on Mh amounts to ~ 130 GeV for Mt = 175 GeV. For the 
two values of tgP introduced above, the Higgs masses Mh, Mh and Mfj± are presented in 
Figs. [pD-d as a function of the pseudoscalar mass M4 for vanishing mixing parameters. 
The dependence on the mixing parameters fi, At, A^ is rather weak and the effects on the 
masses are limited by a few GeV [0. 

The MSSM couplings of Table |I] are shown in Fig. ^ as functions of the pseudoscalar 
mass Ma for two values of tg/3 = 1.5 and 30 and vanishing mixing parameters. The mixing 
effects are weak and thus phenomenologically unimportant. For large values of tg(3 the 
Yukawa couplings to (up) down-type quarks are (suppressed) enhanced and vice versa. 
Moreover, it can be inferred from Fig. ^ that the couplings of the light scalar Higgs particle 
approach the SM values for large pseudoscalar masses, i.e. in the decoupling regime. Thus 
it will be difficult to distinguish the light scalar MSSM Higgs boson from the SM Higgs 
particle, in the region where all Higgs particles except the light scalar one are very heavy. 

1.3 Organization of the Paper 

In this work we will review and update all Higgs decay widths and branching ratios as 
well as all relevant Higgs boson production cross sections at the LHC within the SM and 



MSSM. Previous reviews can be found in Refs. P3| , |3^ . However, this work contains 



substantial improvements due to our use of new results. Moreover, we will use recent 
parametrizations of parton densities for the production cross sections at the LHC. 

This paper is organized as follows. In Section 2 we will review the decay rates and 
production processes of the SM Higgs particle at the LHC. Section 3 will present the cor- 
responding decay rates and production cross sections for the Higgs bosons of the minimal 
supersymmetric extension. A summary will be given in Section 4. 

2 Standard Model 
2.1 Decay Modes 

The strength of the Higgs-boson interaction with SM particles grows with their masses. 
Thus the Higgs boson predominantly couples to the heaviest particles of the SM, i.e. W, Z 
gauge bosons, top and bottom quark. The decays into these particles will be dominant, 
if they are kinematically allowed. All decay modes discussed in this section are obtained 
by means of the FORTRAN program HDECAY |3|, HQ. 



^The program can be obtained from http://wwwcn.cern.ch/'^mspira/ 



8 



50 



50 



20 - 



10 



5 - 



tgP = 30 



hdd: -sina/cosP 




50 100 200 



500 1000 



20 
10 
5 



tgP = 30 
Hdd: cosa/cosP 



tgp=1.5 



50 100 200 



500 1000 




0.1 

0.05 
0.03 



huu: cosa/sinP 



: / tgp=3o 



M, [GeV] 



_1_L 



_1_ 



_l L_l_ 



0.5 



0.2 



0.1 - 



0.05 



-1 1 1 — I — I I I I 



tgp=1.5 



Huu: -sina/sinP 



tgP = 30 



M, [GeV] 



50 100 200 



500 



0.03 
1000 50 



100 200 



500 1000 




Figure 2: The coupling parameters of the neutral MSSM Higgs bosons as a function of 
the pseudoscalar mass Ma for two values of tgj3 = 1.5, 30 and vanishing mixing. They are 
defined in Table ^ 
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2.1.1 Lepton and heavy quark pair decays of the SM Higgs particle 

In lowest order the leptonic decay width of the SM Higgs boson is given by [|10| 



with P = {1 — Amf / M]jY^'^ being the velocity of the leptons. The branching ratio of decays 
into T leptons amounts to about 10% in the intermediate mass range. Muonic decays can 
reach a level of a few 10~^, and all other leptonic decay modes are phenomenologically 
unimportant. 




Figure 3: Typical diagrams contributing to H ^ QQ at lowest order and one-, two- and 
three-loop QCD. 



For large Higgs masses the particle width for decays to b, c quarks [directly coupling 
to the SM Higgs particle] is given up to three-loop QCD corrections [typical diagrams are 
depicted in Fig. R[ by the well-known expression 



T[H QQ] = mliMH) [Aqcd + A,] (7) 
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with 



^QCD 



1 + 5.67 "-^^^^ + (35.94 - 1.36iV^) ( 



TT 



TT 



+ (164.14 - 25.77Nf + 0.259iV|; 



H 



TT 



A, 



2 r 



TT 



1.57 --log 



M 



H 



M? 



in the MS renormahzation scheme; the running quark mass and the QCD couphng are 
defined at the scale of the Higgs mass, absorbing in this way large mass logarithms. The 
quark masses can be neglected in general, except for heavy quark decays in the threshold 
region. The QCD corrections in this case are given, in terms of the quark pole mass Mq, 
by 11 



T[H^QQ 



4v/27r 



1 + — ^A^ 

3 TT 



(8) 



where (3 = {1 — AMq/ Mfj)^/"^ denotes the velocity of the heavy quarks Q. To leading 
order, the QCD correction factor reads as p8[ 



1 (3 + 34/5^-13/5^) log + -1^(7/3^-1) 



1-/5 8/32 



with 



(9) 



4Li2 



l + /5^ 



+ 211, 



1 + /5 



— 3 log log 



/5 °l + /5 



-2 log ^3"^ log /5 



-3/5 log. 



-4/51og/5. 



[Li2 denotes the Spence function, Li2(x) = — Jq dyy ^ log(l— y).] Recently the full massive 
two-loop corrections of 0{Npa1) have been computed; they are part of the full massive 



two-loop result [41 



The relation between the perturbative pole mass Mq of the heavy quarks and the MS 
mass mQ^Mo) at the scale of the pole mass can be expressed as |^ 



m, 



Q 



{Mq) 



M, 



Q 



3 TT ^ \ Tl 



2 ' 



(10) 



where the numerical values of the NNLO coefficients are given by Kt ~ 10.9, ~ 12.4 
and Kc ~ 13.4. Since the relation between the pole mass Mc of the charm quark and the 
MS mass mdMc) evaluated at the pole mass is badly convergent the running quark 
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masses mq^MQ) have to be adopted as starting points. [They have been extracted directly 
from QCD sum rules evaluated in a consistent 0{as) expansion [^].] In the following we 
will denote the pole mass corresponding to the full NNLO relation in eq. ( ]T0| ) by Mq^ and 
the pole mass corresponding to the NLO relation [omitting the contributions of Kq\ by 
Mq*^ according to Ref . [^ . Typical values of the different mass definitions are presented 
in Table 0. It is apparent that the NNLO correction to the charm pole mass is comparable 
to the NLO contribution starting from the MS mass. 



Q 




MS' 




mQ(100 GeV) 


c 


1.23 GeV 


1.42 GeV 


1.64 GeV 


0.62 GeV 


h 


4.23 GeV 


4.62 GeV 


4.87 GeV 


2.92 GeV 


t 


167.4 GeV 


175.0 GeV 


177.1 GeV 


175.1 GeV 



Table 2: Quark mass values for the MS mass and the two different definitions of the 
pole masses. The strong coupling has been chosen as as{Mz) = 0.118, and the bottom 
and charm mass values are taken from Ref. U^] . The last column shows the values of the 
running MS masses at a typical scale fi = 100 GeV. 



The evolution from Mq upwards to a renormalization scale fi can be expressed as 

c[asifi)/TT] , 



mqi^i) = mQ^Mq) 



c[as{Mq)/7:] 



with the coefficient function 44, 45 



c{x) 
c(x) 
c(x) 



^ a;^ [1 + 0.895X + 1.371 + 1.952 



[1 + 1.014X + 1.389 x^ + 1.091 x'^ 



2 

25 

12 

— x\ [1 + 1.175X + 1.501 x^ + 0.1725x^1 

6 / 

7 \t 

-x\ [1 + 1.398X+ 1.793x^- 0.6834x^1 



for Ms < II <Mc 
for Mc < fi <Mb 

for Mb < fi < Mt 

for Mt < jJi. 



Ill 



For the charm quark mass the evolution is determined by eq. (0) up to the scale /i = Mb, 
while for scales above the bottom mass the evolution must be restarted at Mq = Mb- The 
values of the running b, c masses at the scale n = 100 GeV, characteristic of the relevant 
Higgs masses, are typically 35% (60%) smaller than the bottom (charm) pole masses M^*^ 
(M^*^) as can be inferred from the last column in Table ^. Thus the QCD corrections 
turn out to be large in the large Higgs mass regime reducing the lowest order expression 
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[in terms of the quark pole masses] by about 50% (75%) for bottom (charm) quarks. The 
QCD corrections are moderate in the threshold regions apart from a Coulomb singularity 
at threshold, which however is regularized by the finite heavy quark decay width in the 
case of the top quark. 

In the threshold region mass effects are important so that the preferred expression for 
the heavy quark decay width is given by eq. (§). Far above the threshold the massless 
0{al) result of eq. (|^ fixes the most improved result for this decay mode. The transition 
between the two regions is performed by a linear interpolation as can be inferred from 
Fig. 1^, thus yielding an optimized description of the mass effects in the threshold region 
and the renormalization group improved large Higgs mass regime. 
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Figure 4: Interpolation between the full massive NLO expression (dashed line) for the bb 
decay width of the Standard Higgs boson and the renormalization group improved NNNLO 
result (dotted line). The interpolated curve is presented by the full line. 



Electroweak corrections to heavy quark and lepton decays are well under control 
47[| . In the intermediate mass range they can be approximated by 



elw 




log 



Mi; 



kfM^ + M, 



w 



-5 + 



log 



Mi 



6vj — a 



f 



(12) 

with Vf = 2I^f — Aefs'^ and aj = 2/3^. L^j denotes the third component of the electroweak 
isospin, Cf the electric charge of the fermion / and sw = sm9w the Weinberg angle; a 
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denotes the QED coupling, Mj the top quark mass and M^y the W boson mass. The 
large logarithm log Mfj/M'j can be absorbed in the running fermion mass analogous to 
the QCD corrections. The coefficient kj is equal to 7 for decays into leptons and light 
quarks; for h quarks it is reduced to 1 due to additional contributions involving top quarks 
inside the vertex corrections. Recently the two- and three-loop QCD corrections to the 
kf terms have been computed by means of low-energy theorems The results imply 
the replacements 

as{Mt)\ 



ki 



kfX {I 



1 /3 ^ 



for / ^ 6 



A;fe X 1 - 4 (1 + C2 



as{Mt 



(13) 



The three-loop QCD corrections to the kf term can be found in ||5^. The electroweak 
corrections are small in the intermediate mass range and can thus be neglected, but we 
have included them in the analysis. However, for large Higgs masses the electroweak 
corrections may be important due to the enhanced self-coupling of the Higgs bosons. In 



the large Higgs mass regime the leading contributions can be expressed as |^ 
with the coupling constant 



2.12A 



32.66P} 



A 



(14) 



(15) 



16V27r2 ■ 

For Higgs masses of about 1 TeV these corrections enhance the partial decay widths by 
about 2%. 

In the case of ti decays of the Standard Higgs boson, below-threshold decays H 
ti* — * tbW~ into off-shell top quarks may be sizeable. Thus we have included them below 



the tt threshold. Their Dalitz plot density reads as |52 

dr 



-{H ^tt* Wth) 



327r'^ y( + -ftKt 



(16) 



dxidx2 

with the reduced energies Xi^2 = "^Et^/MH and the scaling variables yi^2 = 1 — a^i,2, 
Ki = M'l /Mfj and the reduced decay widths of the virtual particles 7^ = Vl/Mfj. The 
squared amplitude can be written as 

To = ?/i(l - yi - ?/2 + - S/tt) + 2Kvy (2/12/2 - Kvy - 2Ki2/i + 4/ttKH/) 

-f^tVm + /«<(! - 4/tt)(2?/i + Kw + nt) . (17) 

The differential decay width in eq. ([T6|) has to be integrated over the xi,X2 region, which 
is bounded by 

2(1 - Xi - X2 + Kt + Kb - Kw ) + X1X2 ^ ^ 



/Xf - 4:Kt^Jxi - AKb 

The transition from below to above the threshold is provided by a smooth cubic interpo- 
lation. Below-threshold decays yield a tt branching ratio far below the per cent level for 
Higgs masses Mh ^ 2Mt. 
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2.1.2 Higgs decay into gluons 




Figure 5: Diagrams contributing to H ^ gg at lowest order. 



The decay of the Higgs boson into gluons is mediated by heavy quark loops in the Standard 
Model, see Fig. ^; at lowest order the partial decay width |]TD|, Q ^ ^ is given by 



Tlo [H gg] 



36 V2 7r3 



(19) 



with the form factor 



-r[l + (l-r)/(r)] 



arcsm 



log j== - i-n 

1 - y/1 - T 



T > 1 



T < 1 



(20) 



The parameter tq = AMq/Mjj is defined by the pole mass Mq of the heavy loop quark 



H 




rrrsir^ g 



H 



fTRTO^ g 



-mnr^ g H 



LQ_QJL> g 





ULfiiL^ g 



Figure 6: Typical diagrams contributing to the QCD corrections to H ^ gg. 



Q. For large quark masses the form factor approaches unity. QCD radiative corrections 
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are built up by the exchange of virtual gluons, gluon radiation from the quark triangle 
and the splitting of a gluon into two gluons or a quark-ant iquark pair, see Fig. |^. If all 
quarks u, - ■ ■ ,b are treated as massless at the renormalization scale ~ Mh ~ 100 GeV, 
the radiative corrections can be expressed as p3|-p5[] 



r^- [H gg (g), qqg] = T^o [«f ^H^^,)] |l + 

7 

4 6 ^ ^ 



TT 



(21) 



with Np = 5 light quark flavors. The full massive result can be found in [^S]. The 
radiative corrections are plotted in Fig. ^ against the Higgs boson mass. They turn out 
to be very large: the decay width is shifted by about 60-70% upwards in the intermediate 
mass range. The dashed line shows the approximated QCD corrections defined by taking 
the coefficient i?^^ in the limit of a heavy loop quark Q as presented in eq. (^). It can be 
inferred from the figure that the approximation is valid for the partial gluonic decay width 
within about 10% for the whole relevant Higgs mass range up to 1 TeV. The reason for the 
suppressed quark mass dependence of the relative QCD corrections is the dominance of 
soft and coUinear gluon contributions, which do not resolve the Higgs coupling to gluons 
and are thus leading to a simple rescaling factor. 

Recently the three-loop QCD corrections to the gluonic decay width have been evalu- 
ated in the limit of a heavy top quark They contribute a further amount of (9(20%) 
relative to the lowest order result and thus increase the full NLO expression by 0(10%). 
The reduced size of these corrections signals a significant stabilization of the perturbative 
result and thus a reliable theoretical prediction. 

The QCD corrections in the heavy quark limit can also be obtained by means of a 
low-energy theorem |Ty, 0. The starting point is that, for vanishing Higgs momentum 



Ph —>■ 0, the entire interaction of the Higgs particle with W, Z bosons and fermions can 
be generated by the substitution 



Mi 



M- X 



1 + ^^ 

V . 



f,W,Z) 



(22) 



where the Higgs field H acts as a constant complex number. At higher orders this substi- 
tution has to be expressed in terms of bare parameters |]53|, |58|. Thus there is a relation 



between a bare matrix element with and without an external scalar Higgs boson [X de- 
notes an arbitrary particle configuration]: 



lim M{XH) 



1 d 

-moj—MiX). 



(23) 



In most of the practical cases the external Higgs particle is defined as being on-shell, 
so that pjj = Mfj and the mathematical limit of vanishing Higgs momentum coincides 
with the limit of small Higgs masses. In order to calculate the Higgs coupling to two 
gluons one starts from the heavy quark Q contribution to the bare gluon self-energy 
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Figure 7: The size of the QCD correction factor for H gg, defined as T = T lo{^ + ■ 
The full line corresponds to the full massive result, while the dashed line shows the heavy 
top quark limit. The top and bottom masses have been chosen as Mt = 175 GeV, Mb = 
5 GeV and the NLO strong coupling constant is normalized as as{Mz) = 0.118. 



M.{gg). The differentiation with respect to the bare quark mass mg can be replaced by 
the differentiation by the renormahzed MS quark mass TfiQimQ). In this way a finite 
contribution from the quark anomalous mass dimension 7m(as) arises: 



mo 



d rnQirriQ) 



d 



dmo 1 + 7m (as) dmQ{mQ) 



(24) 



The remaining mass differentiation of the gluon self-energy results in the heavy quark 
contribution (3Q{as) to the QCD (3 function at vanishing momentum transfer and to an 
additional contribution of the anomalous dimension of the gluon field operators, which 
can be expressed in terms of the QCD (3 function |5^. The final matrix element can be 
converted into the effective Lagrangian |^3|, ^ ^ |5^, |9l| 



a^*^) [mt(mt) 



4 

and al 



(25) 



a 



^^^fmt(mj)]. The strong coupling as of the effective 



with as 

theory includes only Np = 5 flavors. The effective Lagrangian of eq. 



is valid for the 
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limiting case M'jj <^ 4Mq. The anomalous mass dimension is given by |60 



.a. 



7m(as) = 2— + 



/lOl 
VT2" 



-\Nf + 1\ 



(26) 



Up to NLO the heavy quark contribution to the QCD (3 function coincides with the 
corresponding part of the MS result. But at NNLO an additional piece arises from 
a threshold correction due to a mismatch between the MS scheme and the result for 



vanishing momentum transfer ||6T|-|64 



n-MS, , 11 27 - 2Np a' 

Po 



72 



6 



at 



19 a, 7387- 325A^F ^ ^ 

1 H H ■ ■ 

4 71 288 



TT 



(27) 



The strong coupling constant a, of eq. (|25| ) includes 6 flavors, and its scale is set by 
the top quark mass mt{mt). In order to decouple the top quark from the couplings in 
the effective Lagrangian, the six-flavor coupling a^f^ has to be replaced by the five-flavor 
expression a^^\ They are related by ||63|-|65|pl 



11 / af^[mt{mt)] 



72 



TT 



(28) 



Finally the perturbative expansion of the effective Lagrangian can be cast into the form 



a 



(5) 



H 



V 




• /II /3i\ af(M,) 

U /3o/ TT 
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Po 
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(29) 



where we have introduced the top quark pole mass Mt. The coefficients of the QCD (3 
function in eq. (|29D are given by 



Pi 
P2 



33 - 2Nf 
12 

153 - 19Np 



24 



1 



128 



2857 



5033 



-Nf + 



325 



-Ni 



(30) 



9 27 

^It should be noted that eq. ( [2^ ) differs from the result of Ref. However, the difference can be 
traced back to the Abelian part of the matching relation, which has been extracted by the author from 
the analogous expression for the photon self-energy p6[ . 
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[The four-loop contribution has also been obtained recently [§2].] Np denotes the number 
of light quark flavors and will be identified with 5. For the calculation of the heavy 
quark limit given in eq. ( pi]) the effective coupling has to be inserted into the blobs of 
the effective diagrams shown in Fig. |^. After evaluating these effective massless one-loop 
contributions the result coincides with the explicit calculation of the two-loop corrections 
in the heavy quark limit of eq. (|2lD at NLO. 




Figure 8: Typical effective diagrams contributing to the QCD corrections to H gg in 
the heavy quark limit. 



Using the discussed low-energy theorem, the electroweak corrections of 0{GfM^) to 
the gluonic decay width, which are mediated by virtual top quarks, can be obtained easily. 
For this purpose the leading top mass corrections to the gluon self-energy have to be com- 
puted. The result has to be differentiated by the bare top mass and the renormalization 
will be carried out afterwards. The final result leads to a simple rescaling of the lowest 
order decay width [B7 



TiH ^ gg) = Tlo{H ^ gg) 



1 + 



GpMl 



(31) 



They enhance the gluonic decay width by about 0.3% and are thus negligible. 




Figure 9: Typical diagram contributing to H ^ QQg- 



The final states H —>■ bbg and cog are also generated through processes in which the b, c 
quarks directly couple to the Higgs boson, see Fig. ^. Gluon splitting g ^ bb in H ^ gg 
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increases the inclusive decay probabilities T{H — > bb + . . .) etc. Since b quarks, and 
eventually c quarks, can in principle be tagged experimentally, it is physically meaningful 
to consider the particle width of Higgs decays to gluon and light u, d, s quark final jets 
separately. If one naively subtracts the final state gluon splitting contributions for b and c 
quarks and keeps the quark masses finite to regulate the emerging mass singularities, one 
ends up with large logarithms of the b, c quark masses [in the limit of heavy loop quark 
masses Mq] 

7 1 



6E- 



b,c 



Ml 



log + log 



M2 



(32) 



which have to be added to the hb and cc decay widths [the finite part emerges from the non- 
singular phase-space integrations] . On the other hand the KLN theorem ensures that 
all final-state mass singularities of the real corrections cancel against a corresponding part 
of the virtual corrections involving the same particle. Thus the mass-singular logarithms 
log Mfj/M^^ in eq. (|32|) have to cancel against the corresponding heavy quark loops in the 
external gluons, i.e. the sum of the cuts 1, 2, 3 in Fig. 10 has to be finite for small quark 
masses Mq. [The blobs at the Hgg vertices in Fig. |1^ represent the effective couplings 
in the heavy top quark limit. In the general massive case they have to be replaced by 
the top and bottom triangle loop^.] Thus in order to resum these large final-state mass 
logarithms in the gluonic decay width, the heavy quarks Q = b,c have to be decoupled 
from the running strong coupling constant, which has to be defined with three light flavors, 
if b, c quark final states are subtracted. 




Figure 10: Cut diagrams, involving heavy quark Q loops, contributing to the imaginary 
part of the Higgs self-energy at the two-loop level. 



^It should be noted that the bottom quark triangle loop develops a logarithmic behaviour cx 
M^/Mfj X log^Mfj/M^, which arises from the integration region of the loop momentum, where the 
b quark, exchanged between the two gluons, becomes nearly on-shell. These mass logarithms do not 
correspond to final-state mass singularities in pure QCD and are thus not required to cancel by the KLN 
theorem. 
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H 



af\MH) = af^Mn) x 1 + 



(3)/ 



af\M, 



(33) 



Expressed in terms of three light flavors, the gluonic decay width is free of explicit mass 
singularities in the bottom and charm quark masses. The resummed contributions of 
6, c quark final states are given by the difference of the gluonic widths [eq. (|2l|)] for the 
corresponding number of flavors Np PB[, 



5V[H ^cc+ ..] 
5T\H ^hh + ...] 



p5_p4 



(34) 



in the limit Mjj > 



ML. 



In this way large mass logarithms log M^/M^^ in the remaining 



gluonic decay mode are absorbed into the strong coupling by changing the number of 
active flavors according to the number of contributing flavors in the final states. It should 
be noted that by virtue of eqs. (^) the large logarithms are implicitly contained in the 
strong couplings for different numbers of active flavors. The subtracted parts may be 
added to the partial decay widths into c and h quarks. In a'f'^Mz) the contribution of 
the h quark is subtracted and in af\Mz) the contributions of both the h and c quarks 
are. The values for a'f\Mz) are typically 5% smaller and those of af\Mz) about 15% 
smaller than af\Mz)^ see Table |^. 



afKMz) 


afKMz) 


afKMz) 


0.112 


0.107 


0.101 


0.118 


0.113 


0.105 


0.124 


0.118 


0.110 



Table 3: Strong coupling constants as{Mz) for different numbers of flavors contributing 
to the scale dependence. In a^^^ the b quark contribution is subtracted and in a^^^ the b 
and c quark contributions are. 



2.1.3 Higgs decay to photon pairs 



The decay of the Higgs boson to photons is mediated by W and heavy fermion loops in 
the Standard Model, see Fig. Ill]; the partial decay width [ETf can be cast into the form 



T[H-. 77] 



Gf Ml 
128^2 



Y.Ncfe)A^{rf)+A^M 
f 



(35) 
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Figure 11: Typical diagrams contributing to H —>■ ^yy at lowest order. 



with the form factors 

Afir) = 2r[l + (l-r)/(r)] 
<(r) = _[2 + 3r + 3r(2-r)/(r)] 

and the function /(r) defined in eq. (|20|). The parameters Xj = AMf /Mfj {i = /, W) are 
defined by the corresponding masses of the heavy loop particles. For large loop masses 
the form factors approach constant values: 

Af ^ ioiMl<t:AMl 

^ -7 forM^<4M2, (36) 

The W loop provides the dominant contribution in the intermediate Higgs mass range, 
and the fermion loops interfere destructively. Only far above the thresholds, for Higgs 
masses Mh ~ 600 GeV, does the top quark loop become competitive, nearly cancelling 
the W loop contribution. 



fWWW 7 



H 



t,b 



7 



Figure 12: Typical diagram contributing to the QCD corrections to H 77. 



In the past the two-loop QCD corrections to the quark loops have been calculated 
They are built up by virtual gluon exchange inside the quark triangle [see 
Owing to charge conjugation invariance and color conservation, radiation of a 



Fig. p 
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single gluon is not possible. Hence the QCD corrections simply rescale the lowest order 
quark amplitude by a factor that only depends on the ratios of the Higgs and quark 



masses 



CniTQ) 



-1 



for Mfj < AMI 



(37) 



According to the low-energy theorem discussed before, the NLO QCD corrections in the 
heavy quark limit can be obtained from the effective Lagrangian [p3l p8| 



C 



a 



4 l + 7n 



H 



[a. 



flU 



(Si 



where P^/a = 2(a/7r)[l + a^/vr + ■ ■ ■] denotes the heavy quark Q contribution to the 
QED jS function and 7m (^s) the anomalous mass dimension given in eq. (P^). The NLO 



expansion of the effective Lagrangian reads as [^3|, p8 



2 _^ pfiiy p 2^ 



H 

V 



l-- + 0{ai 

TT 



(39) 



which agrees with the C-value of eq. (^) in the heavy quark limit. 

The QCD corrections for finite Higgs and quark masses are presented in Fig. as a 
function of the Higgs mass. In order to improve the perturbative behaviour of the quark 
loop contributions they should be expressed preferably in terms of the running quark 
masses mQ{MH/2), which are normalized to the pole masses Mq via 



mgifiQ = Mq 



Mr 



(40) 



their scale is identified with fig = within the photonic decay mode. These def- 

initions imply a proper definition of the QQ thresholds Mh = '2'Mq, without artificial 
displacements due to finite shifts between the pole and running quark masses, as is the 
case for the running MS masses. It can be inferred from Fig. |T3| that the residual QCD 
corrections are moderate, of (9(10%), apart from a broad region around M^ ~ 600 GeV, 
where the W loop nearly cancels the top quark contributions in the lowest order decay 
width. Consequently the relative QCD corrections are only artificially enhanced, and the 
perturbative expansion is reliable in this mass region, too. Since the QCD corrections are 
small in the intermediate mass range, where the photonic decay mode is important, they 
are neglected in this analysis. Recently the three-loop QCD corrections to the effective 
Lagrangian of eq. ( pOj) have been calculated |7^. They lead to a further contribution of 
a few per mille. 

The electroweak corrections of 0{GfM^) have been evaluated recently. This part of 
the correction arises from all diagrams, which contain a top quark coupling to a Higgs 
particle or would-be Goldstone boson. The final expression results in a rescaling factor 
to the top quark loop amplitude, given by [|71| 



1 - ^ ( 4:eteb + 5 



14 ,\GfMI 
3^' J 8V27T^ 



(41) 
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Figure 13: The size of the QCD correction factor for H 77, defined hyV = TLoi^ + S) ■ 
The top and bottom masses have been chosen as Mt = 175 GeV, Mb = 5 GeV and the 
strong coupling constant has been normalized to as{Mz) = 0.118 at NLO. The quark 
masses are replaced by their running masses at the scale fiq = Mh/2. 



where Ct^t are the electric charges of the top and bottom quarks. The effect is an en- 
hancement of the photonic decay width by less than 1%, so that these corrections are 
negligible. 

In the large Higgs mass regime the leading electroweak corrections to the W loop 
have been computed by means of the equivalence theorem [0, 173] . This ensures that for 



large Higgs masses the dominant contributions arise from longitudinal would-be Gold- 
stone interactions, whereas the contributions of the transverse W and Z components are 
suppressed. The final result decreases the W form factor by a finite amount [|73| , 



1-3.027 



8V27r2 



for > AMI 



w 



(42) 



These electroweak corrections are only sizeable in the region around Mh ~ 600 GeV, 
where the lowest order decay width develops a minimum due to the strong cancellation 
of the W and t loops and for very large Higgs masses Mh ~ 1 TeV. Since the photonic 
branching ratio is only important in the intermediate mass range, where it reaches values 
of a few 10~^, the electroweak corrections are neglected in the present analysis. 
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2.1.4 Higgs decay to photon and Z boson 
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Figure 14: Typical diagrams contributing to H ^ Z'-f at lowest order. 



The decay of the Higgs boson to a photon and a Z boson is mediated by W and heavy 
fermion loops, see Fig. [T^; the partial decay width can be obtained as |74| 



T[H^ Zj] 



M 



2\3 



64 TT^ 



Ml, 



Y^A'l{rfAj) + A^{Tw.Xw) 
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with the form factors 



2N, 
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ef{hf ~ 2e/ sin^ 6 
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A^{t,\) = cos^H-{4(3-tan2^vF)/2(r,A) 
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1 + - ) tan^ Ow 
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The functions Ji, I2 are given by 
/i(r,A) : 



tX t^X^ 
+ 



2(r-A) 2(r-A)2 
h{r,X) = -^^^^[fi^)-fW] 

where the function g{T) can be expressed as 



[fir) - /(A)] + 



T^X 



(r - A)^ 



[9ir)-g{X)] 
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\Jt — \ arcsin —=z 

IT 



log 
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1 - v^r^ 



r > 1 

r < 1 



(43) 



(44) 



(45) 



and the function /(r) is defined in eq. (pOl). The parameters Xj = AMf/M^ and Aj = 
4Mf /M| (i = /, VT) are defined in terms of the corresponding masses of the heavy loop 
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Figure 15: Typical diagram contributing to the QCD corrections to H ^ Z7. 



particles. Due to charge conjugation invariance, only the vectorial Z coupling contributes 
to the fermion loop so that problems with the axial 75 coupling do not arise. The W loop 
dominates in the intermediate Higgs mass range, and the heavy fermion loops interfere 
destructively. 

The two-loop QCD corrections to the top quark loops have been calculated ||75[ in 
complete analogy to the photonic case. They are generated by virtual gluon exchange 
inside the quark triangle [see Fig. |T5[. Due to charge conjugation invariance and color 
conservation, radiation of a single gluon is not possible. Hence the QCD corrections can 
simply be expressed as a rescaling of the lowest order amplitude by a factor that only 
depends on the ratios Xj and Aj (i = /, W), defined above: 



DHirQ^Q) ^ -1 forM|<M^<4M2. (46) 



In the limit Mz the quark amplitude approaches the corresponding form factor of 
the photonic decay mode [modulo couplings], which has been discussed before. Hence 
the QCD correction in the heavy quark limit for small Z masses has to coincide with 
the heavy quark limit of the photonic decay mode of eq. (|37D . The QCD corrections for 



finite Higgs, Z and quark masses are presented in ||75[ as a function of the Higgs mass. 
They amount to less than 0.3% in the intermediate mass range, where this decay mode 
is relevant, and can thus be neglected. 

2.1.5 Intermediate gauge boson decays 

Above the WW and ZZ decay thresholds, the partial decay widths into pairs of massive 
gauge bosons (V = W, Z) at lowest order [see Fig. |l6l are given by 



T(H ^ VV) = 5vj^m - 4x + 12x') , (47) 



with X = M^/Mfj, (3 = VI -4x and 6v = 2 (1) for V = W (Z). 
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Figure 16: Diagram contributing to H —>■ VV [V = W, Z]. 



The electroweak corrections have been computed in |7g] at the one-loop leveL They 
are small and amount to less than about 5% in the intermediate mass range. Furthermore 
the QCD corrections to the leading top mass corrections of OiGpM^) have been calculated 
up to three loops. They rescale the WW, ZZ decay widths by ^% 



T{H ZZ) 
T{H WW) 



Tlo{H^ZZ)^1 



Xt 



Tlo{H ^WW)\l-Xt 



5-(15-2C2)- 

TT _ 

5-(9-2C2)- 



(48) 
(49) 



with Xt = GpMf /{S\/2'k'^). The three-loop corrections can be found in [^. Since the 
electroweak corrections are small in the intermediate mass regime, they are neglected in 
the analysis. For large Higgs masses, higher order corrections due to the self-couplings of 
the Higgs particles are relevant. They are given by [I79| 



V{H VV) = Tlo{H VV) {l + 2.80A + 62.03^} 



(50) 



with the coupling constant A defined in eq. (p!5|). They start to be sizeable for ^ 400 
GeV and increase the decay width by about 20% at Higgs masses of the order of ~ 1 TeV. 

Below threshold the decays into off-shell gauge particles are important. The partial 
decay widths into single off-shell gauge bosons can be obtained in analytic form ^ 



r(H VV*] 
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'Mf 



(51) 



with 5'^ = l^5'z = 7/12 - 10 sin^ Ow/^ + 40 sin^ ^vf/27 and 
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(52) 



For Higgs masses slightly larger than the corresponding gauge boson mass the decay 
widths into pairs of off-shell gauge bosons play a significant role. Their contribution can 
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be cast into the form 



81 



r(H V*V*\ = 1 f " dQjMyVy riMn-Q.r dQjMyVy 

^ ^ ^ TT^io {Ql-MlY + M^Tlh (Qi-M2)2 + M2r^ ° 

(53) 

with Qi) Q\ being the squared invariant masses of the virtual gauge bosons, My and Vy 
their masses and total decay widths; Fq is given by 



T, = 5y^^^\{QlQlMj,] 



Mi 



(54) 



with the phase-space factor A(x, y\ z) = {1 — x / z — y/zY — ^xyj z^ . The branching ratios 
of double off-shell decays reach the per cent level for Higgs masses above about 100 (110) 
GeV for off-shell W{Z) boson pairs. They are therefore included in the analysis. 

2.1.6 Three-body decay modes 

The branching ratios of three-body decay modes may reach the per cent level for large 
Higgs masses p2[. The decays H — *• W~^W~'y,ti'y{g) are already contained in the QED 



(QCD) corrections to the corresponding decay widths H — > W~^W~ ,ti. However, the 
decay modes H — > W~^W~ Z,tiZ are not contained in the electroweak corrections to the 
WW, tt decay widths. Their branching ratios can reach values of up to about 10~^ for 
Higgs masses M^ ~ 1 TeV. As they do not exceed the per cent level, they are neglected 
in the present analysis. The analytical expressions are rather involved and can be found 



in 82 . 



2.1.7 Total decay width and branching ratios 

In Fig. Il^ the total decay width and branching ratios of the Standard Model Higgs boson 



are shown as a function of the Higgs mass. For Higgs masses below ~ 140 GeV, where 
the total width amounts to less than 10 MeV, the dominant decay mode is the bb channel 
with a branching ratio up to ~ 85%. The remaining 10-20% are supplemented by the 
T~^T~,cc and gg decay modes, the branching ratios of which amount to 6.6%, 4.6% and 
6% respectively, for Mh = 120 GeV [the bb branching ratio is about 78% for this Higgs 
mass]. The 77 [Zj) branching ratio turns out to be sizeable only for Higgs masses 
80 (120) GeV < Mh < 150 (160) GeV, where they exceed the lO'^ level. 

Starting from Mh ~ 140 GeV the WW decay takes over the dominant role joined 
by the ZZ decay mode. Around the WW threshold of 150 GeV < Mh < 180 GeV, 
where the W pair of the dominant WW channel becomes on-shell, the ZZ branching 
ratio drops down to a level of ~ 2% and reaches again a branching ratio ~ 30% above the 
ZZ threshold. Above the ti threshold Mh = 2Mt, the ti decay mode opens up quickly, 
but never exceeds a branching ratio of ~ 20%. This is caused by the fact that the leading 
WW and ZZ decay widths grow with the third power of the Higgs mass [due to the 
longitudinal W, Z components, which are dominating for large Higgs masses] , whereas 
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Figure 17: (a) Total decay width (in GeV) of the SM Higgs boson as a function of its 
mass, (b ) Branching ratios of the dominant decay modes of the SM Higgs particle. All 
relevant higher order corrections are taken into account. 
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the tt decay width increases only with the first power. Consequently the total Higgs 
width grows rapidly at large Higgs masses and reaches a level of ~ 600 GeV at Mh = 1 
TeV, rendering the Higgs width of the same order as its mass. At = 1 TeV the WW 
{ZZ) branching ratio approximately reaches its asymptotic value of 2/3 (1/3). 

2.2 Higgs Boson Production at the LHC 

2.2.1 Gluon fusion: gg ^ H 

The gluon-fusion mechanism ITSf 

pp-^ gg ^ H 

provides the dominant production mechanism of Higgs bosons at the LHC in the entire 
relevant Higgs mass range up to about 1 TeV. As in the case of the gluonic decay mode, 
the gluon coupling to the Higgs boson in the SM is mediated by triangular loops of top 
and bottom quarks, see Fig. |T8[ Since the Yukawa coupling of the Higgs particle to heavy 
quarks grows with the quark mass, thus balancing the decrease of the amplitude, the form 
factor reaches a constant value for large loop quark masses. If the masses of heavier quarks 
beyond the third generation are fully generated by the Higgs mechanism, these particles 
would add the same amount to the form factor as the top quark in the asymptotic heavy 
quark limit. Thus gluon fusion can serve as a counter of the number of heavy quarks, 
the masses of which are generated by the conventional Higgs mechanism. On the other 
hand, if these novel heavy quarks will not be produced directly at the LHC, gluon fusion 
will allow to measure the top quark Yukawa coupling. This, however, requires a precise 
knowledge of the cross section within the SM with three generations of quarks. 




Figure 18: Diagrams contributing to gg ^ H at lowest order. 



The partonic cross section. Fig. |18|, can be expressed by the gluonic width of the Higgs 

CTLoigg ^ H) = 0-0^(1 - z) (55) 



boson at lowest order 15 
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where the scahng variables are defined as 2; = Mjj/s, tq = AMq/M^, and s denotes the 
partonic cm. energy squared. The amphtudes Aq{tq) are presented in eq. (pOl) . 

In the narrow-width approximation the hadronic cross section can be cast into the 
form ||T5[| 

(^Lo{pp ^ H) = aoTH- — (56) 

dTH 



drag i-i dr 

' ''%{x,M')g{T/x,M') (57) 



with 

dr Jt X 

denoting the gluon luminosity at the factorization scale M, and the scaling variable is 
defined, in analogy to the Drell-Yan process, as th = Mfjjs, with s specifying the total 
hadronic cm. energy squared. 




QCD corrections. In the past the two-loop QCD corrections to the gluon-fusion cross 
section. Fig. have been calculated [0, |5^, They consist of virtual corrections 

to the basic gg ^ H process and real corrections due to the associated production of the 
Higgs boson with massless partons. 



99 



Hg and gq — ^ Hq, qq Hg . 



These subprocesses contribute to the Higgs production at 0{al). The virtual corrections 
rescale the lowest-order fusion cross section with a coefficient depending only on the 
ratios of the Higgs and quark masses. Gluon radiation leads to two-parton final states 
with invariant energy s > mjj in the gg, gq and qq channels. The final result for the 
hadronic cross section can be split into five parts [ [55| , 
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with the renormalization scale in as and the factorization scale of the parton densities to 
be fixed properly. The lengthy analytic expressions for arbitrary Higgs boson and quark 
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masses can be found in Refs. [53, 83 1. The quark-loop mass has been identified with the 
pole mass Mq, while the QCD coupling is defined in the MS scheme. We have adopted 
the MS factorization scheme for the NLO parton densities. 

The coefficient C{tq) denotes the finite part of the virtual two-loop corrections. It 
splits into the infrared part vr^, a logarithmic term depending on the renormalization scale 
fj, and a finite quark-mass- dependent piece c{tq), 



Ci^Q) = tt" + c(rQ) + 
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(59) 



The term c{tq) can be reduced analytically to a one-dimensional Feynman-parameter 
integral, which has been evaluated numerically [p3| , p3| , |8^ . In the heavy-quark limit 
tq = AMq/M^j- ^ 1 and in the light-quark limit tq <^ 1, the integrals could be solved 
analytically. 

The finite parts of the hard contributions from gluon radiation in gg scattering, gq 
scattering and qq annihilation depend on the renormalization scale /i and the factorization 
scale M of the parton densities: 
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with z = th/t = Mfj/s] Pgg and Pgq are the standard Altarelli-Parisi splitting functions 
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F, denotes the usual + distribution: F{z)^ = F{z)-5{l-z)!^dz'F{z'). The coefficients 
dgg, dgq Qxid dqg 08.11 be reduced to one-dimensional integrals, which have been evaluated 
numerically Q for arbitrary quark masses. They can be calculated analytically 

in the heavy- and light-quark limits. 

In the heavy-quark limit tq ^ 1 the coefficients c{tq) and dij{z,TQ) reduce to very 



simple expressions EB, 59 



c(^q) ^ y 



dgg{z,TQ) ^-y(l-z)^ 



dgg{z, Tq) 



(l-z)' 



dgq{z, Tq) 



32 
27 



(1-.)= 



(62) 



32 



The corrections of 0{Mfj/M^) in a systematic Taylor expansion have been demonstrated 
to be very small |8^. In fact, the leading term provides an excellent approximation up 
to the quark threshold ~ 2Mq. In the opposite limit where the Higgs mass is much 



larger than the top mass, the analytic result can be found in ||53| 
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Figure 20: K factors of the Q CD- corrected gluon-fusion cross section a{pp —>■ H + X) 
at the LHC with cm. energy y/s = 14 TeV. The dashed lines show the individual contri- 
butions of the four terms of the QCD corrections given in eq. The renormalization 
and factorization scales have been identified with the Higgs mass, = M = Mh and the 
CTEQ4 parton densities have been adopted. 



We define K factors as the ratio 



tot 



(^NLO 



(63) 



The cross sections aNLO in next-to-leading order are normalized to the leading-order 
cross sections (Jlo: convoluted consistently with parton densities and in leading order; 



the NLO and LO strong couplings are chosen from the CTEQ4 parametrizations ||8^ of 
the structure functions, af^°(M^) = 0.116, a^^{Mz) = 0.132. The K factor can be 
decomposed into several characteristic components: K^h accounts for the regularized 
virtual corrections, corresponding to the coefficient C; Kab = g,q,q] for the real 

corrections as defined in eqs. (|6D[). These K factors are presented for LHC energies 
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in Fig. 20 as a function of the Higgs boson mass. Both the renormahzation and the 
factorization scales have been identified with the Higgs mass, fi = M = M^- Apparently 
Kyirt and Kgg are of the same size, of order 50%, while Kgg and Kgg turn out to be quite 
small. Apart from the threshold region Mh ~ 2Mi, Ktot is insensitive to the Higgs mass. 
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Figure 21: Comparison of the exact and approximate NLO cross section a{pp H + 
X) at the LHC with cm. energy ^/s = 14 TeV. The solid line shows the exact cross 
section including the full t, b quark mass dependence and the dashed line corresponds to 
the approximation defined in eq. 1^6^). The renormalization and factorization scales have 
been identified with the Higgs mass, fi = M = Mh and the CTEQ4 parton densities [SljJ 
with NLO strong coupling [as{Mz) = 0.116/ have been adopted. The top mass has been 
chosen as Mt = 175 GeV and the bottom mass as = 5 GeV. 



The corrections are positive and large, increasing the Higgs production cross section at 
the LHC by about 60% to 90%. Comparing the exact numerical results with the analytic 
expressions in the heavy-quark limit, it turns out that these asymptotic K factors provide 
an excellent approximation even for Higgs masses above the top-decay threshold. We 
explicitly define the approximation by 
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where we neglect the b quark contribution in /^^^^^(oo), while the leading order cross 
section includes the full t, b quark mass dependence. The comparison with the full 
massive NLO result is presented in Fig. |21|. The solid line corresponds to the exact cross 
section and the broken line to the approximate one. For Higgs masses below ~ 1 TeV, the 
deviations of the asymptotic approximation from the full NLO result are less than 10%. 

Theoretical uncertainties in the prediction of the Higgs cross section originate from two 
sources, the dependence of the cross section on different parametrizations of the parton 
densities and the unknown NNLO corrections. 
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Figure 22: Higgs production cross section for three different sets of parton densities 
[CTEQ4M, MRS(Rl) and GRV('92)J. 



The uncertainty of the gluon density causes one of the main uncertainties in the 
prediction of the Higgs production cross section. This distribution can only indirectly 
be extracted through order effects from deep-inelastic lepton-nucleon scattering, or 
by means of complicated analyses of final states in lepton-nucleon and hadron-hadron 
scattering. Adopting a representative set of recent parton distributions P^, |SS|, we find 



a variation of about ±10% of the cross section for the entire Higgs mass range. The cross 
section for these different sets of parton densities is presented in Fig. 22 as a function 
of the Higgs mass. The uncertainty will be smaller in the near future, when the deep- 
inelastic electron/positron-nucleon scattering experiments at HERA will have reached the 
anticipated level of accuracy. 
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Figure 23: The renormalization and factorization scale dependence of the Higgs pro- 
duction cross section at lowest and next-to-leading order for two different Higgs masses 
Mh = 150 and 500 GeV. 
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The [unphysical] variation of the cross section with the renormahzation and factoriza- 
tion scales is reduced by including the next-to-leading order corrections. This is shown 
in Fig. ^ for two typical values of the Higgs mass, = 150, 500 GeV. The renormal- 
ization/factorization scale /i = M is varied in units of the Higgs mass, /i = ^Mh for ^ 
between 1/2 and 2. The ratio of the cross sections for ^ = 1/2 and 2 is reduced from 1.62 
in leading order to 1.32 in next-to- leading order for Mfj = 500 GeV. Since, for small Higgs 
masses, the dependence on /i for ^ ~ 1 is already small at the LO level, the improvement 
by the NLO corrections is less significant for a Higgs mass Mh = 150 GeV. However, 
the figures indicate that further improvements are required, because the /i dependence 
of the cross section is still monotonic in the parameter range set by the natural scale 
^ ~ M ~ Mff. The uncertainties due to the scale dependence appear to be less than 
about 15%. 

Soft gluon resummation. Recently soft and coUinear gluon radiation effects for the 
total gluon-fusion cross section have been resummed. The perturbative expansion of the 
resummed result leads to an approximation of the three-loop NNLO corrections of the 
partonic cross section in the heavy top mass limit, which approximates the full massive 
NLO result with a reliable precision [see Fig. |2T[]. Owing to the low-energy theorem 
discussed before [see the gluonic decay mode H — > gg], the unrenormalized partonic cross 
section factorizes, in n = 4 — 2e dimensions, as 



where the coefficient ctq is defined in eq. (j^) with the strong coupling replaced by the 
bare one, aso- The bare correction factor po{z, M]j/ fj,"^, e) arises from the effective diagrams 
analogous to Fig. ^ in higher orders. In the following we shall neglect the contributions 
from gq and qq initial states, which contribute less than ~ 10% to the gluon-fusion cross 
section at NLO. The hadronic cross section can be obtained by convoluting eq. (|65D with 
the bare gluon densities. 
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with the scahng variables z = th /{X1X2) and th = Mfj/s, where s denotes the hadronic 
cm. energy squared. The moments of the hadronic cross section factorize into three 
factors: 

a{N,Ml,fi') = J^TH r^-' a{TH, M],, f,') = gl{N + 1) ^°^(iV, M^, /x^, e) . (69) 



The bare correction factor po may be expanded perturbatively, 
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51, |5|, m, see eq. (|6g): 



The first two [unrenormahzed] coefficients are known from the exphcit NLO calculation 
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where we have absorbed trivial constants into a redefinition of the scale, /x^ — > /x^ exp[7£; — 
log(47r)]. 

The starting point for the resummation is provided by the Sudakov evolution equation 
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(73) 

which follows from the basic factorization theorems for partonic cross sections into soft, 
coUinear and hard parts at the boundaries of the phase space The solution for the 
moments of eq. (^) is given by 
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(74) 



where we have imposed the boundary condition 

Po 



U-f = 0,a,(/i),ej =6{l~z) 



(75) 



which is valid in n dimensions. The bare evolution kernel Wq{z, X,as{fi),e) can be eval- 
uated perturbatively. After renormalizing the strong coupling as and the gluon densities 
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in the MS scheme all singularities cancel, and the finite renormalized correction factor 



reads as 
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with /3o = (33— 2A'"i?)/6. It should be noted that in the last exponential we have kept terms 
of O {\og^ N/N) {i > 1) in the moments of the correction factor, which are not covered 
by the basic factorization theorems near the soft and collinear edges of phase space. On 
the other hand at NLO they turn out to originate from collinear gluon radiation and 
are thus universal, so that they can be included in the resummationQ. In order to define 
the resummed correction factor we have to perform a regularization of the singularity at 
A ~ Aqcd, which is related to an infrared renormalon. Nevertheless, the perturbative 
expansion is well defined. The NLO and NNLO results for fi = M read [91] 
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^ Their inclusion in the Drell-Yan process and deep-inelastic scattering yields the correct coefficients 
of the log"^ N/N terms and those log^ N/N terms, which are related to the strong coupling constant, at 
NNLO, which supports the consistency of their resummation. However, a rigorous proof has not been 
worked out so far. 
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The novel contributions of 0(log* N/N) appear as the non- infrared functions Si{z). They 
are of significant importance for processes at the LHC and therefore have to be included 
to gain a reliable approximation by means of soft gluon resummation. 
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Figure 24: Exact and approximate two- and three-loop correction factor convoluted with 
NLO gluon densities in the heavy top quark limit for the SM Higgs boson. The CTEQ4M 
parton densities have been adopted with as{Mz) =0.116 at NLO. 

The convolution of the correction factor with NLO gluon densities and strong cou- 



pling is presented in Fig. as a function of the Higgs mass at the LHC. The solid line 
corresponds to the exact NLO result and the lower dashed line to the NLO expansion of 
the resummed correction factor. It can be inferred from this figure that the soft gluon 
approximation reproduces the exact result within ~ 5% at NLO. The upper dashed line 
shows the NNLO expansion of the resummed correction factor. From the analogous anal- 
ysis of the Drell-Yan process at NNLO we gain confidence that the NNLO expansion 
of the resummed result reliably approximates the exact NNLO correction [^. Fig. ^ 
demonstrates that the correction factor amounts to about 2-2.3 at NLO and 2.7-3.5 at 
NNLO in the phenomenologically relevant Higgs mass range Mh ^ 1 TeV. However, in 
order to evaluate the size of the QCD corrections, each order of the perturbative expan- 
sion has to be computed with the strong coupling and parton densities of the same order, 
i.e. LO cross section with LO quantities, NLO cross section with NLO quantities and 
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NNLO cross section with NNLO quantities. This consistent K factor amounts to about 



1.5-1.9 at NLO and is thus about 50-60% smaller than the result in Fig. |2^. Therefore 
a reliable prediction of the gluon-fusion cross section at NNLO requires the convolution 
with NNLO parton densities, which are not yet available. Thus it is impossible to predict 
the Higgs production cross section with NNLO accuracy until NNLO structure functions 
will be accessible. 

The scale dependence of the gluon-fusion cross section [neglecting gq and qq initial 
states] is presented in Fig. ^ as a function of the scale in units of the Higgs mass, 
^ = ix/Mh- All orders of the perturbative expansion are evaluated with NLO parton 
densities and strong coupling, so that the LO and NNLO curves do not correspond to 
physically consistent values. The dotted line represents the LO and the lower full line 
the exact NLO scale dependence. The two dashed curves correspond to the NLO and 
the NNLO expansions of the resummed cross section. The upper solid line shows the full 
NNLO scale dependence, which has been obtained from the exact NLO result by means of 

This curve has been identified with the approximate 
I supports the vahdity of the resummed expression 



renormalization group methods |]91 
NNLO expansion at ^ 



1. Fig. _ 

within a reasonable accuracy for physically relevant scale choices 1/2 < ,^ < 2. Moreover, 
the upper solid line clearly indicates that the NNLO scale dependence develops a broad 
maximum around the natural scale /i ~ Mh for large Higgs masses and thus a significant 
theoretical improvement. 



Electroweak corrections. The electroweak corrections to the gluon-fusion cross sec- 
tion have been computed in two different limits. The leading top mass corrections of 
0{GfM^) coincide with the corrections to the gluonic decay mode of eq. (^) and are 
thus small For large Higgs masses the electroweak corrections of OlGpMfj) have 

been evaluated by means of the equivalence theorem They enhance the cross section 
by about 10-20% for large Higgs masses. 



2.2.2 Vector-boson fusion: 



qq 



qqV*V* — > qqH 



The second important Higgs production channel at the LHC is the vector-boson-fusion 



mechanism [see Fig. 26], which will be competitive with the dominant gluon-fusion mech- 
anism for large Higgs masses ~ 1 TeV [jl6|, |l^. For intermediate Higgs masses the 



vector-boson-fusion cross section is about one order of magnitude smaller than the gluon 
one. The leading order partonic vector-boson-fusion cross section |1^ can be cast into 
the form [V = W, Z]: 
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Figure 25: ^ca/e dependence of the Higgs production cross section as a function of the 
common renormalization and factorization scale in units of the Higgs mass for two values 
of Mh — 150,500 GeV. All orders of the cross section are evaluated with NLO parton 
densities [CTEQ4M] and strong coupling constant [as{Mz) — 0.116/. 
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Figure 26: Diagram contributing to qq qqV*V* qqH at lowest order. 



F2{xi,M^)F,{x2,M^) 
Piqi 

F2{xuM^)F2{x2,M^) 

(Pigi)(P2g2) 



pp (Pigi)(P2gi) 

-T l-T 2 5 

qi 

(Pigi)(P2g2)(gi92) 

2 2 
91^2 



Piq2 



Piqi ' 

-^qiq2 



qi 

{P2q2){Piq2 
ql 



Fs{x,,M^)F,{x2,M^ 
2(Pigi)(P2g2) 
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where dPS^ denotes the three-particle phase space of the final-state particles, Mp the 
proton mass, Pi,2 the proton momenta and qi^2 the momenta of the virtual vector bosons 
V*. The functions Pj(x, M^) {i = 1,2,3) are the usual structure functions from deep- 
inelastic scattering processes at the factorization scale M: 

F,{x,M') = Y.i< + ^lM^^M') + q{x,M')] 
g 

F2{x,M^) = 2xJ2{vl + alMx,M') + q{x,M')] 



F-Ax.M' 



-q{x,M^) + q{x,M^)] 



where Vg (a^) are the (axial) vector couplings of the quarks q to the vector bosons V: 
Vq = —Ug = -\/2 for V" = and Vg = 2/3^ — 4eqSin^ 9w, ciq = 2/35 for V = Z. I^q is the 
third weak isospin component and Cg the electric charge of the quark q. 

In the past the QCD corrections have been calculated within the structure function 
approach [|r^. Since, at lowest order, the proton remnants are color singlets, no color will 
be exchanged between the first and the second incoming (outgoing) quark line and hence 
the QCD corrections only consist of the well-known corrections to the structure functions 
Pj(x, M^) {i = 1,2,3). The final result for the QCD-corrected cross section leads to the 
replacements 

P,(a;,M2) ^ P,(x,M2) + AP,(x,M2,Q2) (i = 1,2,3) 
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where z = x/y and the functions P^g, P^^ denote the well-known Altarelli-Parisi splitting 
functions, which are given by |85 

4 , . . 3 
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The physical scale Q is given by Q"^ = —qf for x = Xi {i = 1,2). These expressions have to 
be inserted in eq. (^) and the full result expanded up to NLO. The typical renormalization 
and factorization scales are fixed by the vector-boson momentum transfer fi = M = Q. 
The K factor, defined as i^' = (Jnlo/<^lo, is presented in Fig. |^ as a function of the 
Higgs mass. The size of the QCD corrections amounts to about 8-10% and is thus small 



2.2.3 Higgs-strahlung: qq V* VH 



The Higgs-strahlung mechanism qq ^ V* ^ VH {V = W, Z) [see Fig. |28| may be 
important in the intermediate Higgs mass range due to the possibility to tag the associated 
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Figure 27: K factor of the QCD corrections to VV H as a function of the SM Higgs 
mass. The CTEQ4M parton densities have been adopted, and the running strong coupling 
constant has been normalized to as{Mz) = 0.116 at NLO. 



vector boson. Its cross section is about one to two orders of magnitude smaller than the 
gluon-fusion cross section for Higgs masses < 200 GeV. The lowest-order partonic 



cross section can be expressed as [IS 



^Loiqq VH) 
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(1 - M2 /Q2)2 



^6) 



where X{x,y;z) = (1 — x/z — y/z)"^ — Axy/z^ denotes the usual two-body phase-space 
factor and Vq (ag) are the (axial) vector couplings of the quarks q to the vector bosons V, 



which have been defined after eq. (plf). The partonic cm. energy squared s coincides at 
lowest order with the invariant mass = My^ of the Higgs-vector-boson pair squared, 
s = Q^. The hadronic cross section can be obtained from convoluting eq. (|86D with the 
corresponding (anti) quark densities of the protons: 



c^Loipp qq-^ VH) 



-'to a 



da'i 
~d^ 



cTLoiQ'^ = rs) , 



(87) 
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Figure 28: Diagram contributing to qq ^ V* ^ VH at lowest order. 



with To = {Mh + MvY / s and s the total hadronic cm. energy squared. 

The QCD corrections are identical to the corresponding corrections to the Drell-Yan 
process. They modify the lowest order cross section in the following way pO[| : 

a{pp VH) = (TlO + ^(^qq + ^CFqg 

Acr„5 = / dry— — / dz aLo[Q = rzs) ujgq{z) 

n Jto g dr Jto/t 

Aag, = ^fdrY^^f dz aLoiQ' = rzs) uqg{z) (88) 

^ ■^O q^q Jto/t 

with the coefficient functions 

LUqqiz) = -P,,(z)log— + -|2[C2-2]5(l-^)+4Pi(^)-2(l + ^)log(l-^)| 

1 / \ 1 

^..(^) = --Pq,{z)\og\^-^^—-^^j+-{l + 6z-7z'}, (89) 

where M denotes the factorization and fi the renormalization scale. The natural scale 
of the process is given by the invariant mass of the Higgs-vector-boson pair in the final 
state, fi = M = Q. The K factors, defined as K = <7nlo/(^lo, are shown in Fig. |2^ as 
a function of the Higgs mass. The size of the QCD corrections is of about 25-40% and 



they are thus of moderate magnitude |20 



2.2.4 Higgs bremsstrahlung off top quarks 

In the intermediate mass range the cross section of the associated production of the Higgs 
boson with a tt pair can reach values similar to those of the Higgs-strahlung process. It 
may thus provide an additional possibility to find a Higgs boson with mass Mh ^130 
GeV by tagging the additional ti pair and the rare photonic decay mode — > 77 [0. 
This process takes place through gluon-gluon and qq initial states at lowest order [see 
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Figure 29: K factor of the QCD corrections to V* —>■ HV as a function of the SM Higgs 
mass. The CTEQ4M parton densities have been adopted, and the running strong coupling 
constant has been normalized to as{Mz) = 0.116 at NLO. The solid line corresponds to 
W bremsstrahlung and the dashed to Z bremsstrahlung. 



Fig. The result for the lowest order cross section is too involved to be presented here. 
We have recalculated the cross section and found numerical agreement with Refs. |l^, ^ . 

At the LHC the gluon-gluon channel dominates due to the enhanced gluon structure 
function analogous to the leading Higgs production mechanism via gluon fusion. The 
QCD corrections to the Hti production are still unknown. They require the evaluation 
of several one-loop five-point functions for the virtual corrections and real contributions 
involving four particles in the final state, where three of them [H, t, t\ are massive. 

2.2.5 Cross sections for Higgs boson production at the LHC 

The results for the cross sections of the various Higgs production mechanisms at the LHC 
are presented in Fig. |3l|, which is an update of Ref. as a function of the Higgs mass. 
The total cm. energy has been chosen as ^/s = 14 TeV, the CTEQ4M parton densities 
have been adopted with as{Mz) = 0.116, and the top and bottom masses have been set to 
Mt = 175 GeV and Mf, = 5 GeV. For the cross section of Htt and Hbb we have used the 
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Figure 30: Typical diagrams contributing to qq/gg —>■ Htt at lowest order. 



leading order CTEQ4L parton densities due to the fact that the NLO QCD corrections 
are unknown. Thus the consistent evaluation of this cross section requires LO parton 
densities and strong coupling. The latter is normalized as as{Mz) = 0.132 at lowest 
order. The gluon-fusion cross section provides the dominant production cross section for 
the entire Higgs mass region up to Mh ~ 1 TeV. Only for Higgs masses Mh ^ 800 GeV 
the yy-fusion mechanism qq —>■ qqH becomes competitive and deviates from the gluon- 
fusion cross section by less than a factor 2 for Mh ^ 800 GeV. In the intermediate mass 
range the gluon-fusion cross section is at least one order of magnitude larger than all other 
Higgs production mechanisms. At the lower end of the Higgs mass range Mh ^ 100 GeV 
the associated production channels of H + V, H + tt yield sizeable cross sections of about 
one order below the gluon-fusion process and can thus allow for an additional possibility 
to find the Higgs particle. 

The search for the standard Higgs boson will be different within three major mass 
ranges, the lower mass range Mh ^ 140 GeV and the higher one, 140 GeV < Mh ^ 800 
GeV, and the very high mass region Mh ^ 800 GeV. 

Mh < 140 GeV 

In the lower mass range the standard Higgs particle dominantly decays into bb pairs. 
Because of the overwhelming QCD background the signal will be very difficult to extract. 
Only excellent b tagging, which may be provided by excellent /x-vertex detectors, might 
allow a sufficient rejection of the QCD background |^|, although this task seems to be 



very difficult ||2^. The associated production of the Higgs boson with a tt pair or a W 
boson may increase the significance of the H ^ bb decay due to the additional isolated 
leptons from t and W decays, but the rates will be considerably smaller than single Higgs 
production via gluon fusion |l^, |19 



Studies for the detection of the H t^t^ decay mode have also been performed. 
Again because of the overwhelming backgrounds from tt and Drell-Yan t~^t~ pair pro- 
duction, this possibility has been found to be hopeless for the Standard Model Higgs 



boson [94|. The branching ratio into off-shell Z pairs H ^ Z*Z* is too small to allow for 
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Figure 31: Higgs production cross sections at the LHC [y/s = 14 TeV] for the various 
production mechanisms as a function of the Higgs mass. The full QCD-corrected results 
for the gluon fusion gg —>■ H , vector boson fusion qq — > VVqq —>■ Hqq, vector boson 
bremsstrahlung qq V* HV and associated production gg, qq —>■ Hti, Hbb are shown. 
The QCD corrections to the last process are unknown and thus not included. 



a detection of four-lepton final states p5 



The only promising channel for the detection of the Higgs boson with masses Mh ^ 140 
GeV is provided by the rare if — >■ 77 decay mode [l^ with a branching ratio of O{10~^). 
For an LHC luminosity of / £ = 10^pb~^ the cross section times branching ratio for pp — » 
77) +X yields C(0.5-l x 10^) events in the mass range 80 GeV < Mh < 140 GeV. 
In order to reject the large backgrounds from the 77 continuum production and the two- 
photon decay mode of the neutral pions tt" — ^ 77, the detection of the rare photonic 
decay mode requires excellent energy and geometric resolution of the photon detectors 
[0]. Moreover, a necessary rejection factor of 10^ for jets faking photons in the detector 
seems to be feasible |]l4l- Thus the LHC studies conclude that the rare photonic decay 
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mode will be a possibility to find the standard Higgs particle in the lower mass range. 



140 GeV <Mh< 800 GeV 

Above the ZZ threshold, on-shell H —>■ ZZ — ^ 4Z^ decays of the Higgs particle provide a 
very clean signature with small SM backgrounds |^ . The two pairs of electrons or muons 



of this 'gold-plated' decay channel carry invariant masses equal to the Z boson mass, thus 
allowing for very stringent cuts against background processes. Below the ZZ threshold, 
off-shell H — s> ZZ* — ^ 4/^ decays, where one of the Z bosons is on-shell, yield clean 
signatures with rather small SM backgrounds [^. However, in the mass range 155 GeV 
^ Mh ^ 180 GeV, where the ZZ branching ratio drops down to values of about 2%, 
the number of events at the LHC allows for a discovery of the Higgs boson only, if the 



maximal luminosity will be reached |jT^. On the other hand the dominant WW decay 
mode of the Higgs boson leads to Z+/~z/z/ final states with strong spin correlations of the 
visible charged lepton pair. A recent analysis has shown that the Higgs particle can easily 



be detected within a few days in this mass range ||23|| . 
Mh > 800 GeV 

For large Higgs masses the total Higgs decay width exceeds 100 GeV and reaches a value 
of about 600 GeV for Mh = 1 TeV. Thus the Higgs resonance peaks in the 4-lepton final 
states become broad and, owing to the decreasing number of events with growing Higgs 
mass, the 'gold-plated' signal H ZZ 41^ will no longer be visible. In order to extend 
the Higgs search to masses beyond 1 TeV, the decay modes H —>■ ZZ, WW —>■ l^l^vv will 
be the only possible signatures. The present status of the studies is not fully conclusive, 
but promising |T^ . 



Fig. ^ shows the expected signal significance at the LHC as a function of the SM 
Higgs mass after using the full experimental data samples of both experiments, ATLAS 
and CMS. It is apparent that after reaching the full integrated luminosity the SM Higgs 



signal may be extracted in the whole relevant mass region 14 



3 Minimal Supersymmetric Extension of the Stan- 
dard Model 

The couplings of the MSSM Higgs bosons to MSSM particles grow with the MSSM particle 
masses, if these are generated by the Higgs mechanism. Thus the MSSM Higgs bosons 
predominantly couple to heavy quarks and gauge bosons. However, for large values of tgj3 
the couplings to down-type quarks are enhanced, so that the coupling to bottom quarks 
may be much larger than to top quarks. Moreover, the Higgs boson interaction with the 
intermediate gauge bosons is always reduced with respect to the SM. The decays into 
heavy particles will be dominant, if they are kinematically allowed. The analysis includes 
the complete radiative corrections to the MSSM Higgs sector due to top/bottom quark 
and squark loops within the effective potential approach, as discussed in the introduction. 
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Figure 32: Expected significance of the SM Higgs boson search at the LHC as a function 
of the Higgs boson mass after reaching the anticipated integrated luminosity J C = 10^pb~^ 
and combining the experimental data of both LHC experiments, ATLAS and CMS. Pro- 
duced from Refs. - courtesy of F. Cianotti. 



51 



Next-to-leading order QCD corrections and the full mixing in the stop and sbottom sectors 
are incorporated. The corresponding formulae are based on the works of Ref. [0. As 



for the SM case, the decay widths and branching ratios of the MSSM Higgs bosons are 
evaluated by means of the FORTRAN program HDECAY [^. 



3.1 Decay Modes 

3.1.1 Decays into lepton and heavy quark pairs 

At lowest order the leptonic decay width of neutral MSSM Higgs bosonQ decays is given 
by [|10|, g 

r[$ - /^n = ^(,f ) (90) 



where gf denotes the corresponding MSSM coupling, presented in Table ||, /5 = (1 — 
4mf/M|)^/^ the velocity of the final-state leptons and p = 3 (1) the exponent for scalar 
(pseudoscalar) Higgs particles. The r pair decays play a significant role, with a branching 
ratio of up to about 10%. Muon decays can develop branching ratios of a few lO"'^. All 
other leptonic decay modes are phenomenologically irrelevant. 

The analogous expression for the leptonic decays of the charged Higgs reads as 



1= — m,tg p 



mj 



(91) 



The decay mode into r+z/,- reaches branching ratios of about 100% below the tb threshold 
and the muonic one ranges at a few 10~^. All other leptonic decay channels of the charged 
Higgs bosons are unimportant. 

For large Higgs masses [M$ ^ Mq] the QCD-corrected decay widths of the MSSM 
Higgs particles into quarks can be obtained from evaluating the analogous diagrams as 
presented in Fig. |, where the Standard Model Higgs particle H has to be substituted by 
the corresponding MSSM Higgs boson $ ||38|-^: 



r[$ ^ QQ] 



—-^mQ{M^){gQ) 



A 



QCD 



(92) 



Neglecting regular quark mass effects, the QCD corrections Aqcd are presented in eq. (|^) 
and the top quark induced contributions read as [HO 



A 



h/H 



9t 



h/H 



'as{M, 



h/H) 



h/H 

9q 



TT 



1.57- -loj 
3 ^ 



Mi 



1 
9 



+ ^ log' 



h/H) 



A^ = ^ 

9q 



TT 



3.83 



2 ^UMa) 



Ml 



'In the following we denote the different types of neutral Higgs particles by (f> 



K/H 



h, H, A. 
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Analogous to the Standard Model case the large logarithmic contributions of the QCD 
corrections are absorbed in the running MS quark mass mqlM^) at the scale of the cor- 
responding Higgs mass M$. In the large Higgs mass regimes the QCD corrections reduce 
the bb (cc) decay widths by about 50 (75)% due to the large logarithmic contributions. 

The heavy quark decay width of the charged Higgs boson reads, in the large Higgs 
mass regime Mh± ^ Mjj + M^), as | 



97] 



IV? 



UD\ 



TT 



-.A\2 



A 



QCD 



(93) 



[Eq. ( p3D is valid if either the first or the second term is dominant.] The relative couplings 
Qq have been collected in Table |I| and the coefficient Vud denotes the CKM matrix element 
of the transition oi D to U quarks. The QCD correction factor Aqcd is given in eq. (0), 
where large logarithmic terms are again absorbed in the running MS masses rnu^D{MH±) 
at the scale of the charged Higgs mass Mh±- In the large Higgs mass regimes, the QCD 
corrections reduce the cb and cs decay widths by about 50-75%, because of the large 
logarithmic contributions. 

In the threshold regions mass effects play a significant role. The partial decay widths 
of the neutral Higgs bosons ^ = h,H and A into heavy quark pairs, in terms of the quark 



pole mass Mq, can be cast into the form 



r[<i> ^ QQ] 



TT 



1 + l^A* 

3 TT 



(94) 



where /5 = (1 — 4Mq/M|)^/^ denotes the velocity of the final-state quarks and p = 3 (1) 
the exponent for scalar (pseudoscalar) Higgs bosons. To next-to-leading order, the QCD 
correction factor is given by eq. (^ for the scalar Higgs particles h, H, while for the 
CP-odd Higgs boson A they read correspondingly as 



A^ 



(95) 



with the function A{(3) defined after eq. (|^). The QCD corrections in the tt threshold 
region are moderate, apart from a Coulomb singularity, which is regularized by taking 
into account the finite top quark decay width. 

The partial decay width of the charged Higgs particles into heavy quarks may be 
written as p7| 



T[H- 



UD 



I VUDI 



4V27T 




Ml 



1 + 



4 a. 



A 



UD 



AMuMd^JmjJj^ 



3 TT 

1 + -^A 



(96) 



4 a, 



3 TT 



UD 



where /Xj = Mf /M'jj±, and A = (1 — /it/ — /^d)^ — 'ifJ'UfJ^D denotes the usual two-body 
phase-space function; the quark masses M^/xi are the pole masses. The QCD factors 
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with the scahng variables Xi = 2/ij/[l — /ij — /i^ + A^/^] and the generic function 
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The transition from the threshold region, involving mass effects, to the renormalizat ion- 
group- improved large Higgs mass regime is provided by a smooth linear interpolation 
analogous to the SM case in all heavy quark decay modes. 

The full MSSM electroweak and SUSY-QCD corrections to the fermionic decay modes 
have been computed |]^. They turn out to be moderate, less than about 10%. Only for 
large values of tg/5 > 10 do the gluino corrections reach values of 20 to 50%, if the relevant 
squark masses are less than ~ 300 GeV. The electroweak and SUSY-QCD corrections are 
neglected in this analysis. 

Below the tt threshold, heavy neutral Higgs boson decays into off-shell top quarks are 
sizeable, thus modifying the profile of these Higgs particles significantly in this region. 
The dominant below-threshold contributions can be obtained from the SM expression 
eq. ([T|) 



dT 



[H W ^ Wth) = {g, 



H\2 



dV 



SM 



tt* 



Wth). 



(9J 



dxidx2 dxidx2 
The corresponding dominant below-threshold contributions of the pseudoscalar Higgs par- 



ticle are given by ||5/ 
dV 



dxidx2 



(A ^ tt* Wtb) 



327r- 



Tn 



yi + iti^t 



(99) 



with the reduced energies Xi^2 = '^Et^/MA, the scaling variables yi^2 = 1 — 2:1^25 = 
Mf /M\ and the reduced decay widths of the virtual particles 7i = V'j/M\. The squared 
amplitude may be written as |5^ 



To = yl{l - yi - y2 + i^w - i^t) + 2Kw{yiy2 - k^w) - '^t{yiy2 - 2?/i - kw - i^t) ■ (100) 
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The differential decay widtlis of eqs. (pHD, (|99[) liave be integrated over tlie Xi,X2 region, 
bounded by eq. (|I8|) . In tliese formulae W and charged Higgs boson exchange contribu- 
tions are neglected, because they are suppressed with respect to the off-shell top quark 
contribution to Wtb final states. However, for the sake of completeness they are included 
in the analysis. Their explicit expressions can be found in [^]. The transition from be- 
low to above the threshold is provided by a smooth cubic interpolation. Below-threshold 
decays yield a tt branching ratio at the per cent level already for heavy scalar and pseu- 
doscalar Higgs masses Mh^a ~ 300 GeV. 

Below the tb threshold off-shell decays t*b bbW^ are important. For Mh± < 



Mt + Mb — Tf their expression can be cast into the form ||52| 



Y{H+^t*b^Wbb) = ^^\^L mh± l^iAKwKt + 3k, - Akw) log 



+ (3/s:^ - - + 1) log ~ ^ (101) 

{3k^ - KtKw - 2KtK^r + Ak^,) + 4 - -Kw 



with the scaling variables /tj = Mf/M'jj± {i = t,W). The b mass has been neglected 
m eq. ([TOlD , but it has been taken into account in the present analysis by performing 
a numerical integration of the corresponding Dalitz plot density, given in [Q. The off- 
shell branching ratio can reach the per cent level for charged Higgs masses above about 
100 GeV for small tg(3, which is significantly below the tb threshold Mh± ~ 180 GeV. 

3.1.2 Gluonic decay modes 



h,H,A 




Q h,H 



Q 



Figure 33: Typical diagrams contributing to ^ —>■ gg at lowest order. 



Since the b quark couplings to the Higgs bosons may be strongly enhanced for large tg/3 
and the t quark couplings suppressed in the MSSM [see Fig. |^, b loops can contribute 
significantly to the gg coupling so that the approximation Mq ^ Mfj can in general no 
longer be applied. The leading order width for h,H gg is generated by quark and 
squark loops, the latter contributing significantly for squark masses below about 400 GeV 
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99]. The contributing diagrams are depicted in Fig. 33. The partial decay widths are 
given by @, |5|, |99| 
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for « ~ 



and the MSSM 



= Q,Q). The function /(r) is defined in eq. ( 
couphngs gn^ can be found in Table |1|. The squark couplings g~^^ are summarized in 
Table The amplitudes approach constant values in the limit of large loop particle 
masses: 



1 
1 

4 



for 
for 



«4M2 
<4M^. 



The squark loop contributions are significant for squark masses < 400 GeV and neg- 
ligible above This can be inferred from Fig. where the ratio of the gluonic decay 
width with and without the squark contributions is shown as a function of the squark mass 
Mg for two values of tgP = 1.5, 30. The QCD corrections to the squark contribution are 
only known in the heavy squark mass limit. The relative QCD corrections are presented 
in Fig. ^ as a function of the corresponding Higgs mass for two representative values 
of tg/3 = 1.5,30. The solid lines include the top and bottom quark as well as squark 
contributions [for = 200 GeV] and the dashed lines only the quark contributions. 
The comparison of the solid and dashed curves implies that the squark loop contributions 
cause a small effect on the relative QCD corrections, so that a reasonable approximation 
within about 10% to the gluonic decay width can be obtained by multiplying the full 
lowest order expression with the relative QCD corrections including only quark loops. 

In complete analogy to the quark contributions the heavy squark loop correction can 
be obtained by means of the extension of the previously described low-energy theorem to 
scalar squark particles [^. The effective NLO Lagrangian for the squark part is given, 
according to eq. (|25|), by 
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Figure 34: Ratio of the QCD-corrected decay width r{h — >■ gg) with and without squark 
loops for two values of tg(3 = 1.5,30 as a function of the common squark mass M^. 
The pseudoscalar has been identified with Ma = 100 GeV. The secondary axes show the 
corresponding values of the light scalar Higgs mass. 



where Pgids) = «s/(127r)[l + llas/{2n)] denotes the heavy squark contribution to the 
QCD f] function ||100|| and 7m(«s) = 4as/(37r) the anomalous squark mass dimension 
p.01|| . Up to NLO the effective couphng is described by 



487r 



25 a,- 
1 + 

6 TT . 



(105) 



Thus the only difference to the quark loops in the heavy loop mass limit arises in the 
virtual corrections. This leads to the additional last term of eq. ( |103| ). 

It turns out a posteriori that the heavy quark limit M^^^ <^ 4Mq is an excellent 
approximation for the QCD corrections within a maximal deviation of about 10% in the 
parameter ranges where this decay mode is relevant. 



For the pseudoscalar Higgs decays only quark loops are contributing, and we find [53 
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Figure 35: Size of the QCD correction factor for h/H — > gg, defined as T = r^o(l + 
as a function of the corresponding Higgs mass for two values of tgf3 = 1.5,30. The full 
lines include the full mass dependence on the top and bottom masses and, in addition, 
the squark contributions in the heavy-squark limit. The dashed curves correspond to the 
omission of the squark contributions. 
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97 7 
T ~ 6^ 



for Ml < 4M^ 



with tq = AMq/M^. The MSSM couphngs gq can be found in Table |T]. For large quark 
masses the quark amplitude approaches unity. In order to get a consistent result for 
the two-loop QCD corrections, the pseudoscalar 75 coupling has been regularized in the 
't Hooft-Veltman scheme ||102|| , which requires an additional finite renormalization of the 
AQQ vertex [^, |103|| . The relative QCD corrections are presented in Fig. |3^ as a function 
of the pseudoscalar Higgs mass Ma for two values of tg/3 = 1.5, 30. The heavy quark limit 
M\ -C 4Mq provides a reasonable approximation in the MSSM parameter range where 
this decay mode is significant. At the threshold M4 = 2Mt, the QCD corrections develop 
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Table 4: MSSM Higgs couplings to charged Higgs bosons, charginos and sfermions relative 
to SM couplings. Qu and Su {i = 1,2) are related to the mixing angles between the 
charginos xf md xf, see Refs. J^, 



a Coulomb singularity, which will be regularized by including the finite top decay width 



The heavy quark limit can also be obtained by means of a low-energy theorem. The 
starting point is the ABJ anomaly in the divergence of the axial vector current 

2n 



with G^j, = |e^!yo/3(j"^"^ denoting the dual field strength tensor. Since, according to 
the Sutherland- Veltman paradox [106], the matrix element {gg\d^i^^\^) vanishes for zero 



momentum transfer, the matrix element {gg\MQQi'~)^QY^) of the Higgs source can be 
related to the ABJ anomaly in eq. ( |108| ). Thanks to the Adler-Bardeen theorem, the 
ABJ anomaly is not modified by radiative corrections at vanishing momentum transfer 
[[L05|] , so that the effective Lagrangian 

C^fS = g^^G'^^^G%^ (109) 
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Figure 36: Size of the QCD correction factor for A — * gg, defined as T = T^oi^ + 
a function of the pseudoscalar Higgs mass for two values of tgP = 1.5, 30. 



is valid to all orders of perturbation theory. In order to calculate the full QCD corrections 
to the gg decay width, this effective coupling has to be inserted in the effective diagrams 
analogous to those of Fig. |^. The final result agrees with the explicit expansion of the 
two-loop diagrams in terms of the heavy quark mass. 

In analogy to the SM case the bottom and charm final states from gluon splitting may 
be added to the corresponding bb and cc decay modes so that the number of light flavors 
has to be chosen as Np = 3 in the scalar and pseudoscalar decays into gluons |^ . 



3.1.3 Decays into photon pairs 

The decays of the scalar Higgs bosons to photons are mediated by W and heavy fermion 
loops as in the Standard Model and, in addition, by charged Higgs, sfermion and chargino 
loops; the relevant diagrams are shown in Fig. ^ The partial decay widths P, ^ are 
given by 
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Figure 37: Typical diagrams contributing to $ ^ 77 at lowest order. 



(110) 



with the form factors 



4f,(r) = 2r[l + (l-r)/(r)] 



A 



h/H 



V) = -r[l-r/(r)] 



[2 + 3r + 3r(2-r)/(r)] 



where the function /(r) is defined in eq. (pOl). For large loop particle masses the form 
factors approach constant values, 
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3 
1 

3 



for « 4M| 

for « 4MJ^ ~ 

-7 for M^/^ < 4M2. . 



Sfermion loops start to be sizeable for sfermion masses M~ < 300 GeV. For larger sfermion 
masses they are negligible. 

The photonic decay mode of the pseudoscalar Higgs boson is generated by heavy 



charged fermion and chargino loops, see Fig. |3^. The partial decay width reads as P, |53| 

2 



T{A' - 77) 



with the amplitudes 



32v^7r3 
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1111 
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For large loop particle masses the pseudoscalar amplitudes approach unity. 



The parameters = AMf / M| 



/, W, if^, x^, /) are defined by the corresponding 



mass of the heavy loop particle and the MSSM couplings g'^^^^^ _^ are summarized in 
Tables |T] and § ' ' ' 

The QCD corrections to the quark and squark loop contributions have been evaluated. 
For the t, b quark loops they are known for finite quark and Higgs masses while 
in the case of squark loops only the large squark mass limit has been computed so far 
107|| . The QCD corrections rescale the lowest order quark amplitudes |107|| , 



Ch/Hirq) 
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for Ml,H « 4M2 
for Ml < AMI 



(113) 
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h/H\ 



vr 



(114) 



for « 4M 



Q 



The QCD corrections to the 77 decay width are plotted in Fig. ^ for two values of 
tg/3 = 1.5, 30 in the case of heavy charginos and sfermions. They are defined in terms of 
the running quark masses in the same way as the SM photonic decay width. The QCD ra- 
diative corrections are moderate in the intermediate mass range |^3|, ^ , where this decay 
mode will be important, and therefore neglected in the analysis. Owing to the narrow- 
width approximation of the virtual quarks, the QCD corrections to the pseudoscalar decay 
width exhibit a Coulomb singularity at the ti threshold, which is regularized by taking 



into account the finite top quark decay width |104 



The QCD corrections to the quark loops in the heavy quark limit can be obtained 
by means of the low-energy theorems for scalar as well as pseudoscalar Higgs particles, 
which have been discussed before. The result for the scalar Higgs bosons agrees with the 
SM result of eq. (|37D, and the QCD corrections to the pseudoscalar decay mode vanish in 
this limit due to the Adler-Bardeen theorem. In complete analogy to the gluonic decay 
mode, the effective Lagrangian can be derived from the ABJ anomaly and is given to all 
orders of perturbation theory by 



eff 



(115) 



Since there are no effective diagrams generated by light particle interactions that con- 
tribute to the photonic decay width at next-to-leading order, the QCD corrections to the 
pseudoscalar decay width vanish, in agreement with the explicit expansion of the massive 
two-loop result. 

Completely analogous the QCD corrections to the squark loops for the scalar Higgs 
particles in the heavy squark limit can be obtained by the extension of the scalar low- 
energy theorem to the scalar squarks. Their coupling to photons at NLO can be described 
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Figure 38: Size of the QCD correction factor for $ 77, defined as T = T^oi^ + 
as a function of the corresponding Higgs mass for two values of tgP = 1.5,30. The 
renormalization scale of the running quark masses is identified with hq = M$/2. The 
common squark mass has been chosen as Ms — 1 TeV. 
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by the effective Lagrangian ||107 



^eff = 9, 



e2~ 

H Q 



Q 4 l + 7„(a. 



H 



(116) 



where = a/{2n)[l + ias/n] denotes the heavy squaA contribution to the QED P 
function ||100|| and 7m(«s) = 4:as/{37i) the anomalous squark mass dimension [|101|| . Up 
to NLO the effective couphng reads as ||107 



a 



H 

in ■ V 

This correction is small and thus neglected in the present analysis. 
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3.1.4 Decays into Z boson and photon 
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Figure 39: Typical diagrams contributing to $ — > Z7 at lowest order. 



The decays of the scalar Higgs bosons into Z boson and photon are mediated by W 
and heavy fermion loops as in the Standard Model and, in addition, by charged Higgs, 
sfermion and chargino loops; the contributing diagrams are shown in Fig. |3^. The partial 
decay widths read as |5^, |108|| 
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with the form factors ^J^, A^r^ given in eq. (|4^), and 

.h/H, cos26w J , , X 

COS Uw 
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where the function Ji(r, A) is defined after eq. 

The Z7 decay mode of the pseudoscalar Higgs boson is generated by heavy charged 
fermion and chargino loops, see Fig. The partial decay width is given by [p.08|| 



T{A ^ Z7) 
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^1 



Xj Xj Xi Xj 



with the fermion amplitudes 



Af{T,X) 



2N, 



3/ 



2e/ sin^ 9 
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/2(r,A). 



cos 9w 

The contributions of charginos and sfermions involve mixing terms. 



(120) 



:i2ii 



expressions can be found in 



Their analytical 
For large loop particle masses and small Z mass, 



the form factors approach the photonic amplitudes modulo couplings. The parameters 



f, W, H"^, x"^, f) are defined by the corresponding 



7 are 



mass of the heavy loop particle and the non-mixing MSSM couplings g 

summarized in Tables |I| and ^, while the mixing and Z boson couplings gf can be found 
ini. 

The branching ratios of the Z7 decay modes range at a level of up to a few 10 
in the intermediate mass ranges of the Higgs bosons and are thus phenomenologically 
unimportant in the MSSM. 



3.1.5 Decays into intermediate gauge bosons 

The partial widths of the scalar MSSM Higgs bosons into W and Z boson pairs can be 
obtained from the SM Higgs decay widths by rescaling with the corresponding MSSM 
couplings gy^ -I which are listed in Table |l|: 

V{h/H -> yWyW) = {g\/^fV{HsM yWyW) . (122) 

They are strongly decreased by kinematic suppression and reduced MSSM couplings, and 
thus do not play a dominant role as in the SM case. Nevertheless the WW^ ZZ branching 
ratios can reach values of (9(10%) for the heavy scalar Higgs boson H for small tg/3. 
Off-shell WW, ZZ decays can pick up several per cent of the light scalar Higgs decays at 
the upper end of its mass range. The pseudoscalar Higgs particle does not couple to W 
and Z bosons at tree level. 



3.1.6 Decays into Higgs particles 

The heavy scalar Higgs particle can decay into pairs of light scalar as well as pseudoscalar 
Higgs bosons, see Fig. ^ The partial decay widths are given by p 



T{H hh) 



X2 



GpM: 



16V2tiMh\ 



Ml 



(123) 
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Figure 40: Typical diagrams contributing to Higgs decays with Higgs bosons in the final 
state. 



1-4^ (124) 



The self-couplings Xuhh and Xhaa can be derived from the effective Higgs potential [pl[1 . 
The decay mode into pseudoscalar particles is restricted to small regions of the MSSM 
parameter space, where the pseudoscalar mass Ma is small. The decay into light scalar 
bosons is dominant for small ig(3 below the ti threshold. 

The contributions of final states containing off-shell scalar or pseudoscalar Higgs 
bosons may be significant and are thus included in the analysis. Their expressions read 
as 
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H — ^ "'^ ^ f arctan — -^=— = — arctan ] I . (125) 

where = M|/M^. They slightly enhance the regions, where the hh,AA decay modes 
of the heavy scalar Higgs boson H are sizeable. 

Moreover, Higgs bosons can decay into a gauge and a Higgs boson, see Fig. The 
various partial widths can be expressed as 

nn^AZ) = >^HAZ ^^^^^ \/A(M1, M|; M^)A(M1, M^; M|) (126) 

4 



C M' 

TiH^H^W^) = ^HH^w ^^^^^ ^X{Mj,,,M^; M^)A(M^,, M^; M^) (127) 

C I 

T{A^hZ) = XU, ^ I ^X{MlMhMl)X{MlMlMl) (128) 



' Sv^vrM 
GfM, 



A 



4 



T{H+^hW+) = ^Ih+w ^J^^^^ \lKMl M^; M^^)A(M,^ M^,; M^) ,(129) 
where the couplings Xf.f. can be determined from the effective Higgs potential ||3l|. The 



functions X{x, y; z) = {1 — x/z — y/z)"^ — 4:xy/z'^ denote the usual two-body phase-space 
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factors. The branching ratios of these decay modes may be sizeable in specific regions of 
the MSSM parameter space. 

Below-threshold decays into a Higgs particle and an off-shell gauge boson turn out to 
be very important for the heavy Higgs bosons of the MSSM. The individual contributions 



are given by p2| 



T{H AZ*) 


,2 ^,9GlMlMj,^ 


(130) 
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The generic functions Gij can be written as 
1 



Gij = - \ 2(-l + Kj - Ki)^JXij 



TT / kJI — Ki + Kj) — X 

- + arctan ^ ^ ^ 

2 I il-K,).fX. 



using the parameters 



2Ki) log Hi + -(1 - Ki 



5(1 + Ki) - 4:Kj 



1 -|- -|- 2Kj (/?j Kj) , 
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-A,; 



Ml 
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(135) 
(136) 

(137) 



The coefficient 5'^ is defined after eq. (|5T|). Off-shell hZ* decays are important for the 



pseudoscalar Higgs boson for masses above about 130 GeV for small ig(3 The decay 
modes if^ hW*, AW* reach branching ratios of several tens of per cent and lead to 
a significant reduction of the dominant branching ratio into tu final states to a level of 
60-70% for small tg/3 0. 



3.1.7 Total decay widths and branching ratios of non-SUSY particle decays 



Fig. ^ presents the total decay widths and Fig. ^ the branching ratios of the various 
Higgs decay modes into non-SUSY particles, i.e. SM and Higgs particles, as a function of 
the corresponding Higgs masses for two representative values of tgjS = 1.5,30. Since the 
Higgs self-interactions are determined by the gauge couplings, the total decay widths of 
all MSSM Higgs bosons do not exceed about 30 GeV, so that these states will appear as 
rather narrow resonances. The small decay widths are a direct consequence of the absence 
of quadratic divergences in the MSSM Higgs sector and the solution of the hierarchy 
problem. 
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Figure 41: Total decay widths of the MSSM Higgs bosons h,H,A,H^ for non-SUSY 
decay modes as a function of their masses for two values of tg[3 = 1.5, 30 and vanishing 
mixing. The common squark mass has been taken to be Ms — 1 TeV. 
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Figure 42: Branching ratios of the MSSM Higgs bosons h{a), H{b), A{c), H^{d) for non- 
SUSY decay modes as a function of their masses for two values of tgj3 = 1.5,30 and 
vanishing mixing. The common squark mass has been chosen as Ms = 1 TeV. 
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For the light scalar Higgs boson h the bb decays dominate, with a branching ratio 
of up to about 90%, see Fig. |4^a. The bulk of the remaining decay modes is taken by 
r+r" decays, the branching ratio of which ranges at about 8-9%. At the upper bound of 
the light Higgs boson mass all decay modes, as for the intermediate SM Higgs particle, 
are important. Their branching ratios coincide with the SM values for the corresponding 
SM Higgs mass, in accordance with the condition that in the decoupling regime the light 
scalar Higgs particle behaves as the SM Higgs boson. 

Fig. shows that for large tgjS the heavy scalar Higgs boson H predominantly 
decays into bb final states with a branching ratio of about 90%, and to a lesser extent into 
r+r^ pairs with a branching ratio of about 10%. All other decay modes are unimportant 
for large tg/3. In contrast, the heavy scalar Higgs particle exhibits a very rich spectrum 
of decay modes for small tg/5. For tg/? = 1.5 the hh decay mode plays the dominant 
role below the tt threshold with a branching ratio of up to 90%. Only in the vicinity of 
Mh ~ 130 GeV does this decay mode drop down, because the trilinear self-coupling Xuhh 
changes sign and crosses zero. This is the only range where the bb decay channels provides 
the dominant contribution, but it falls off very quickly above and below this Higgs mass. 
Moreover, WW decays are sizeable with a branching ratio of about 10-30% below the tt 
threshold, while the ZZ decays reach values of less than 8%. Above the tt threshold, tt 
decays are overwhelming and their branching ratio amounts to up to 98%. 

From Fig. ^2| c it can be inferred that for large tg/3 the pseudoscalar Higgs particle A 
only decays into bb [BR ~ 90%] and r+r^ pairs [BR ~ 10%]. All other decay channels 
are suppressed and thus unimportant. Contrary to that at small tg/? the bb decay mode 
dominates only below the Zh threshold with a branching ratio ~ 80-90%. The branching 
ratio of t~^t~ decays ranges at about 8-9% in this mass regime. Above the Zh threshold, 
the Zh decay channel plays the dominant role and its branching ratio can reach about 
50% below the ti threshold. It should be noted that already below the Zh threshold 
off-shell Z*h decays are sizeable and thus important. In addition the gg decay channel 
grows rapidly from 2% up to about 20% at the tt threshold. Above this threshold tt 
decays overwhelm with a branching ratio of nearly 100%. 

Fig. ^ d shows that below the tb threshold charged Higgs r^Ur decays provide 

the dominant contribution. Owing to the sizeable below-threshold decays into W*h and 
W*A, the branching ratio of the rur decays does not exceed 70% for small tg/?, but 
amounts to about 99% for large tg/5. Above the tb threshold, if^ —>■ tb is dominant. 
For small tg/5 its branching ratio reaches about 99%, whereas for large tg/3 it does not 
exceed about 80% due to a still sizeable contribution of r^Ur decays. For small tgP a long 
off-shell tail below the tb threshold arises from off-shell t*b decays. Just below the 

tb threshold Wh decays can be dominant for small tgf3 within a very restricted charged 
Higgs mass range. For small charged Higgs masses the off-shell decays into W*h and W*A 
can acquire branching ratios of more than 10% for small tg/3. Below the Wh threshold 
cs and cb decays reach branching ratios of a few per cent. 
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3.1.8 Decays into SUSY particles 



Chargino/neutralino masses and couplings. The chargino/neutralino masses and 
couplings to the MSSM Higgs bosons are fixed by the Higgs mass parameter /i and the 
SU(2) gaugino mass parameter M2. The mass matrix of the charginos is given by ||25 



M2 V2Mw sin P 

\f2M-w cos /5 [i 



(138) 



This can be diagonalized by two mixing matrices f/, V ^ yielding the masses of the physical 
xf 2 states: 



M + 

Al,2 
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w 



I ^ 1/2 

Ty'(M| - /i2)2 + cos2 2(3 + 4M^(M| + /i^ + 2M2/i sin 2/3) ^ (139) 

If either /z or M2 is large, one chargino corresponds to a pure gaugino state and the other 



to a pure higgsino state. The Higgs couplings to charginos ||109|, p.lCI|| can be expressed as 
[k = 1, 2, 3, 4 correspond to H, h, A, H^] 



(140) 



where the coefficients and dk are defined to be 
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— sin [3 
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(141) 



The mass matrix of the four neutralinos depends in addition on the U(l) gaugino mass 
parameter Mi, which is constrained by SUGRA models to be Mi = |tan^vK^2- In the 



bino-wino-higgsino basis, it has the form [p5 



Ml -Mzsmewcos(3 MzsmOwsm/S 

M2 Mz cos 9w cos P —Mz cos 9w sin /5 

—Mz sin ^14/ cos /3 Mz cos cos /3 —fj, 

sin 6'vi/ sin /? — cos6'vi/ sin/3 —/i 

(142) 

which can be diagonalized by a single mixing matrix Z. The final results are too involved 



to be presented here. They can be found in | 109 |. If either or M2 is large, two neutralinos 
are pure gaugino states and the other two pure higgsino states. The Higgs couplings to 
neutralino pairs ||109| , |110|| can be written as [k = 1, 2, 3, 4 correspond to H, h, A, H^] 



Hk XiXj '■ Fijk 
with the coefficients ek,dk defined in eq. (|141|) 



-(Zj2 — tan6'vyZji)(efc^j3 + dkZa) 



(143) 
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The charged Higgs couphngs to chargino-neutrahno pairs are fixed to be ||109 
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UiiZj^ — —j=Ui2{Zj2 + tanOwZji] 



(144) 



Sfermion masses and couplings. The scalar partners fi.R of the left- and right- 
handed fermion components mix with each other. The mass eigenstates /1.2 of the 
sfermions / are related to the current eigenstates fi^R by mixing angles Of, 



fi = fLCOs6f + fn sin 6'/ 
h = -/l sin 6*/ + //J cos 61/ , 



(145) 



which are proportional to the masses of the ordinary fermions. Thus mixing effects are 
only important for the third-generation sfermions t, b, f, the mass matrix of which is given 
by m 

- Ml+Mj Mf{Af-f,rf) 
_Mj{Aj-f,rf) Ml + Mj 

with the parameters rt = r-r = 1/rt = tg/5. The parameters Aj denote the Yukawa mixing 
parameters of the soft supersymmetry breaking part of the Lagrangian. Consequently the 
mixing angles acquire the form 
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and the masses of the squark eigenstates are given by 



M 



The neutral Higgs couplings to sfermions read as ||109 



(147) 
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'IlIr 2 

with the couplings gf listed in Table |[ The charged Higgs couplings to sfermion pairs 



109|| can be expressed as [a, (3 = L,R] 



9lX = -^[9r + M'^g^, 



2 al3-i 



(150) 



with the coefficients g"^ summarized in Table 
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Table 5: Coefficients of the neutral MSSM Higgs couplings to sfermion pairs. 
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Table 6: Coefficients of the charged MSSM Higgs couplings to sfermion pairs. 



Decays into charginos and neutralinos. The decay widths of the MSSM Higgs 
particles [k = 1,2,3,4 correspond to H, h, A, H"^] into neutralino and chargino pairs 
can be cast into the form [p.09| , |110| 
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(151) 



where ?7i,2,4 = +1, = ~1 ^'^^ 



unless the final state consists of two identical 



(Majorana) neutralinos, in which case da = 1; = ±1 stands for the sign of the i'th 
eigenvalue of the neutralino mass matrix, which can be positive or negative. For charginos 
these parameters are always equal to unity. The symbols Xij^k denote the usual two-body 
phase-space functions 
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(152) 



If chargino/neutralino decays are kinematically allowed, which may be the case for the 
heavy MSSM Higgs particles if, A, if^, their branching ratios can reach values up to 



74 



about 100% below the corresponding top quark thresholds. They can thus be dominant, 
jeopardizing the Higgs search at the LHC due to the invisibility of a significant fraction of 
these decay modes ||109|| . A typical example of the total sum of chargino/neutralino decay 
branching ratios is shown in Fig. ^ for the heavy Higgs bosons. Even above the corre- 
sponding top quark thresholds the chargino/neutralino branching ratios will be sizeable. 
For large Higgs masses they reach common values of about 20-80%: in the asymptotic 
regime Mh^. ^ M^, the total sum of decay widths into charginos and neutralinos acquires 
the simple form 



109, 110 
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5Z XiXj 



4v/27r 



1 H — tan^ 6w 
3 



(153) 



for all three Higgs bosons if, A, i/^, which is independent of any MSSM parameter 
[tg/3, /i, At^ft, M2]. Normalized to the total width, which is dominated by ti^bb (tb) decay 
modes for the neutral (charged) Higgs particles the branching ratio of chargino/neutralino 
decays will exceed a level of about 20% even for small and large tg/3. In some part of 
the MSSM parameter space, invisible light scalar Higgs boson decays into the lightest 
neutralino h — > XiXi "wiH be possible and their branching ratio can exceed 50% [|109|, |110| . 



Decays into sleptons and squarks. The sfermionic decay widths of the MSSM Higgs 
bosons Hk [k = 1, 2, 3, 4 correspond to H, h, A, and i,j = 1, 2] can be written as ||109 



^<^'^^-^^'' = vi;s^v^<«M>^- ''''' 

The physical MSSM couplings gf) can be obtained from the couplings presented in 

Jijj 

eqs. (|149|) and ( |150|) by means of the mixing relations in eq. (|145|) . The symbol Xij^k 



denotes the usual two-body phase-space factor of eq. 

In the limit of massless fermions, which is a valid approximation for the first two 
generations, the pseudoscalar Higgs bosons A do not decay into sfermions due to the 
suppression of sfermion mixing by the fermion mass. In the decoupling regime, where the 
Higgs masses Mfj H± are large, the decay widths of the heavy scalar and charged Higgs 
particles into sfermions are proportional to ||109 



T{H, ^ //) oc sin^ 2/3 . (155) 

Thus they are only important for small tg/3 ~ 1. However, they are suppressed by an 
inverse power of the large Higgs masses, rendering unimportant the sfermion decays of 
the first two generations. 

Decay widths into third-generation sfermions [t, b, f] can be much larger, thanks to 
the significantly larger fermion masses. For instance, in the asymptotic regime the heavy 
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Figure 43: Branching ratios of the MSSM Higgs boson H, A, decays into 
charginos/neutralinos and squarks as a function of their masses for tg(3 = 1.5. The mix- 
ing parameters have been chosen as fi = 160 GeV, At = 1.05 TeV, A^ = and the squark 
masses of the first two generations as = 400 GeV. The gaugino mass parameter has 
been set to M2 = 150 GeV. 



scalar Higgs decay into stop pairs of tlie same lielicity is proportional to [ |109| ] 



T{H ^ it) oc 



GfMI_ 



(156) 



which will be enhanced by large coefficients compared to the ffist /second-generation 
squarks for small tg/5. At large tg/3 sbottom decays will be significant. Moreover, for 
large Higgs masses the decay widths of heavy neutral CP-even and CP-odd Higgs parti- 



cles into stop pairs of different helicity will be proportional to ||10£ 



T{H,A-^ it) oc 



GpMl 
Id. 



/i + 



At 



tg/3 



(157) 



and hence will be of the same order of magnitude as standard fermion and chargino/neu- 
tralino decay widths. In summary, if third-generation sfermion decays are kinematically 
allowed, they have to be taken into account. An extreme example for the total branching 
ratios of decays into squarks is depicted in Fig. where they can reach values of ~ 80% 
for the heavy scalar Higgs boson H. 
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Very recently the SUSY-QCD corrections to the stop and sbottom decays of the MSSM 
Higgs bosons have been calculated ||111|| . They reach about 30%, especially in the thresh- 
old regions. They are not included in the present analysis. 



3.2 Neutral Higgs Boson Production at the LHC 
3.2.1 Gluon fusion: 5(51 ^ $ [^ = h,H,A] 

The gluon-fusion mechanism |T3| 

dominates the neutral Higgs boson production at the LHC in the phenomenologically 
relevant Higgs mass ranges for small and moderate values of tgjS. Only for large tg/3 can 
the associated $66 production channel develop a larger cross section due to the enhanced 
Higgs couplings to bottom quarks ||T^, |112|| . Analogous to the gluonic decay modes, the 



gluon coupling to the neutral Higgs bosons in the MSSM is built up by loops involving 



top and bottom quarks as well as squarks, see Fig. 44 




Figure 44: Typical diagram contributing to gg ^ ^ at lowest order. 



The partonic cross sections can be obtained from the gluonic widths of the Higgs 
bosons at lowest order 



(158) 
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where the scaling variables are defined as 2; = M^/s, ti = AM^/M^ {i = Q,Q), and 
s denotes the partonic cm. energy squared. The amplitudes A* ^(tq q) are defined in 
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eqs. ( |T02D , (TO, and the MSSM couplings g%, g% can be found in Tables I and |. In the 
narrow-width approximation the hadronic cross sections are given by 

(TLoipp ^ = (^qT^-; — (159) 

with the gluon luminosity defined in eq. (0) and the scaling variables r$ = M|/s where s 
specifies the total hadronic cm. energy squared. For small tg(3 the top loop contribution 
is dominant, while for large values of tg/5 the bottom quark contribution is strongly 
enhanced. If the squark masses are less than ~ 400 GeV, their contribution is significant. 



and for squark masses beyond ~ 500 GeV they can safely be neglected |99|. This is 
demonstrated in Fig. where the ratio of the cross section with and without the squark 
contribution is presented as a function of the corresponding scalar Higgs mass. In the 
phenomenological mass range the squark loops may enhance the cross section by up to a 
factor 2. 




M, [tgP=1.5] 



' ' ' m5[GeV] 
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Figure 45: Ratio of the Q CD- corrected cross section a{pp —^h + X) with and without 
squark loops as a function of the common squark mass for two values of tg(3 = 1.5, 30, 
and for Ma = 100 GeV. The secondary axes present the corresponding light scalar Higgs 
mass Mfi- The top and bottom masses have been chosen as Mt = 175 GeV, Mf, = 5 GeV, 
and the cross sections are convoluted with CTEQ4M parton densities using as{Mz) = 
0.116 as the normalization of the NLO strong coupling constant. 
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QCD corrections. In the past the two-loop QCD corrections to the gluon-fusion cross 
section were calculated |5^, p.03|| . In complete analogy to the SM case they consist of 
virtual corrections to the basic gg ^ process and real corrections due to the associated 
production of the Higgs bosons with massless partons, 



99 



^g and gq $g, qq <^g 



Thus the contributions to the final result for the hadronic cross section can be classified 
as 

1 + 



a{pp -> $ + X) = (T, 



TT 



firgg 



(160) 



The analytic expressions for arbitrary Higgs boson and quark masses are rather involved 
and can be found in [Q. As in the SM case the (s)quark-loop masses have been identified 
with the pole masses Mq {M~), while the QCD coupling is defined in the MS scheme. We 
have adopted the MS factorization scheme for the NLO parton densities. The axial 75 
coupling has been regularized in the 't Hooft-Veltman scheme ||102|| , which preserves the 
chiral symmetry in the massless quark limit and fulfills the non-renormalization theorem 



of the ABJ anomaly at vanishing momentum transfer |105 



The coefficients C*(rQ, t~) split into the infrared vr^ term, a logarithmic term including 



the renormalization scale /z, and finite (s)quark mass-dependent pieces c'^{tq,t~] 

33 



C*(rQ,r.)=vr^ + c*(rQ,r.) + 
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-log 



Ml 



(161) 



The terms c*(rQ) originating from quark loops have been reduced analytically to one- 

|T0^ . 



dimensional Feynman-parameter integrals, which were evaluated numerically p3t 
The QCD corrections to the squark contributions are only known in the heavy-squark 
limit , which however provides a reasonable approximation to the K factor due to the 
dominance of soft and collinear gluon radiation for heavy particle loops in the gluon-fusion 
process. 

The remaining contributions of eq. ( |160|) can be cast into the form [53, 103 
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with z = Ti^/r = M^/s. Pgg and Pgq are the standard Altarelli-Parisi splitting functions 
defined in eq. (^l]). The coefficients d^g, dgg and d^-^ have been reduced to one-dimensional 
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integrals for the quark loops, which can be evaluated numerically 103|| for arbitrary 



quark masses. They can be calculated analytically in the heavy- and light-quark limits. 

In the heavy-quark limit the quark contributions to the coefficients c'^{tq) and dfj{z, tq) 
reduce to the same expressions as in the SM case of eq. (|62D for the scalar Higgs particles 
h, H. For the pseudoscalar Higgs boson only the coefficient c^{tq) differs from the scalar 
case, 

TQ = AMI/ Ml » 1 ■■ c-^(^q) ^ 6 . 

(163) 

In fact, the leading terms in the heavy-quark limit provide a reliable approximation for 
small tg/5 up to Higgs masses of ~ 1 TeV as can be inferred from Fig. which shows 
the exact pseudoscalar cross sections (solid lines) as a function of the pseudoscalar Higgs 
mass for three values of ig(3 and the approximation obtained by multiplying the full 
massive leading-order cross section with the K factor obtained in the heavy-quark limit. 
A maximal deviation ~ 25% for tg/3 < 5 occurs in the intermediate mass range. The 
squark contribution in the heavy-squark limit coincides with the heavy-quark case apart 
from the virtual piece 



= 4M|/M,V » 1 : c^l-[r^) ^ | . 

(164) 

The QCD corrections to the squark loops have been evaluated for degenerate squark 
masses, so that no mixing effects occur, and for heavy gluinos, such that their contribu- 
tions are suppressed. In this case there are no squark loop effects in pseudoscalar Higgs 
production. 

In the opposite limit, where the Higgs mass is much larger than the quark mass, the 
analytic results coincide with the SM expressions for both the scalar and pseudoscalar 
Higgs particles |53[. This coincidence reflects the restoration of the chiral symmetry in 



the massless quark limit. 

The K factors Ktot = (^nloI^lo are presented for LHC energies in Fig. as a function 
of the corresponding Higgs boson mass. Both the renormalization and the factorization 
scales have been identified with the Higgs masses /i = M = M$. The variation of the 
K factors with the Higgs masses is mainly caused by the MSSM couplings apart from 
the threshold region, where in the pseudoscalar case a Coulomb singularity emerges in 



analogy to the gluonic and photonic decay modes |103| |. The corrections are positive 



and large, increasing the MSSM Higgs production cross sections at the LHC by up to 
about 100%. 

The effect of the squark loops on the scalar Higgs K factors is presented in Fig. 
which shows the K factors of scalar Higgs boson production with and without squark loops 
as a function of the corresponding Higgs mass. It is clearly visible that the squark loops 
hardly change the K factors, making the K factors from the pure quark contributions an 
excellent approximation within maximal deviations of about 10%. 

Theoretical uncertainties in the prediction of the Higgs cross section originate from two 
sources, the dependence of the cross section on different parametrizations of the parton 
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Figure 46: Comparison of the exact and approximate NLO cross section a{pp — > A + X) 
at the LHC with cm. energy s/s = 14 TeV. The solid lines show the exact cross sections 
including the full t, b quark mass dependence and the dashed lines correspond to the heavy- 
quark approximation of the K factor. The renormalization and factorization scales have 
been identified with the Higgs mass, fi = M = M4 and the CTEQ4M parton densities 
with NLO strong coupling [as{Mz) = 0.116/ have been adopted. The top mass has been 
chosen as Mt = 175 GeV, the bottom mass as Mb = 5 GeV and the common squark mass 
as Ms = I TeV. 



densities and the unknown NNLO corrections. For representative sets of recent parton dis- 
tributions we find a variation of about ±10% of the cross section for Higgs masses 
larger than ~ 100 GeV analogous to the SM case. The uncertainty due to the gluon den- 
sity will be smaller in the near future when the deep-inelastic electron/positron-nucleon 
scattering experiments at HERA will have reached the anticipated level of accuracy. 

The [unphysical] variation of the cross sections with the renormalization and factoriza- 
tion scales is reduced by including the NLO corrections. This is shown in Fig. ^ for the 
heavy scalar and pseudoscalar Higgs particles with masses M^ = 500 GeV and Ma = 100 
GeV. The renormalization/factorization scale /x = M is varied in units of the Higgs mass 
/i = ^M$. The remaining uncertainties due to the scale dependence appear to be less 
than about 15%. 
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Figure 47: K factors of the QCD- corrected gluon-fusion cross section a{pp — * $ + X) at 
the LHC with cm. energy y/s = 14 TeV. The dashed lines show the individual contribu- 
tions of the four terms of the QCD corrections given in eq. ( \16Q ). The renormalization and 
factorization scales have been identified with the corresponding Higgs mass, /i = M = M$, 
and the CTEQ4M parton densities have been adopted. 
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Figure 48: K factors of the cross sections a{pp — >• h/H+X) with [solid lines] and without 
[dashed lines] squark loops as a function of the corresponding scalar Higgs mass for two 
values of tg(3 = 1.5,30. The common squark mass has been chosen as = 200 GeV. 
The top and bottom masses have been set to Mt = 175 GeV, Mb = 5 GeV, and the NLO 
cross sections are convoluted with GTEQ^M parton densities using as{Mz) = 0.116 as the 
normalization of the NLO strong coupling constant. The LO cross sections are evaluated 
with GTEQ4L parton densities with the LO strong coupling as{Mz) = 0.132. 



Soft gluon resummation. Recently soft and coUinear gluon radiation effects for the 



total gluon-fusion cross section have been resummed ||9T| . In complete analogy to the SM 
case, the perturbative expansion of the resummed result leads to a quantitative approxi- 
mation of the three-loop NNLO corrections of the partonic cross section in the heavy top 
mass limit, which approximates the full massive NLO result with a reliable precision for 
small and medium values of tg/3 [see Fig. Owing to the low-energy theorems discussed 
before [see the gluonic decay modes $ gg], the unrenormalized partonic cross section 
factorizes in the same way as the SM cross section. The scalar factors k^^^ coincide with 
the SM values of eq. (|66D and the pseudoscalar factor is equal to unity, because of the 
non-renormalization of the ABJ anomaly ||105|| , 

K^ = l. (165) 

The resummation of soft and collinear gluon effects proceeds along the same lines as in 
the SM case. The final results for the scalar correction factors p^^^ are identical to the 
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Figure 49: The renormalization and factorization scale dependence of the Higgs produc- 
tion cross section at lowest and next-to-leading order for two different Higgs bosons H,A 
with masses Mh — 500 GeV and Ma — 100 GeV and two values of tg(5 — 1.5, 30. 
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SM result eq. (|76D, and the pseudoscalar correction factor can be cast into the form |91 

M-i , \ fWMDI ^^^^^ 
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(/i) =p^/^ iV,^,a,(/x) xexp 6- 



TT 



The perturbative expansions at NLO and NNLO |^T| read as [for /i = M] 



P^)h U ^] + 3 {24Pi(2) - 12L,Vo{z) - 48Si{z) 



(167) 



+ (12C2 + 6 + /3oL^)5(l-^)} 
where we use the same notation as in the SM case. 
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Figure 50: Exact and approximate two- and three-loop correction factor convoluted with 
NLO gluon densities in the heavy top quark limit for the pseudoscalar MSSM Higgs boson. 
The CTEQ4M parton densities have been adopted with as{Mz) = 0.116 at NLO. 



The convolution of the scalar correction factors with NLO gluon densities and strong 
coupling coincides with the SM case in Fig. El, while the pseudoscalar case is presented in 
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Fig. 50 as a function of the pseudoscalar Higgs mass at the LHC The sohd hne corresponds 
to the exact NLO result and the lower dashed hne to the NLO expansion of the resummed 
correction factor. It can be inferred from this figure that the soft gluon approximation 
reproduces the exact result within ~ 5% at NLO. The upper dashed hne shows the NNLO 



expansion of the resummed correction factor. Fig. |5y demonstrates that the correction 
factor amounts to about 2.2-2.5 at NLO and 3.2-4.1 at NNLO in the phenomenologically 
relevant Higgs mass range ^ 1 TeV. However, in order to evaluate the size of the 
QCD corrections, each order of the perturbative expansion has to be integrated with the 
strong coupling and parton densities of the same order, i.e. LO cross section with LO 
quantities, NLO cross section with NLO quantities and NNLO cross section with NNLO 
quantities. This consistent K factor amounts to about 1.5-2.0 at NLO and is thus about 



50-60% smaller than the result in Fig. ^ A reliable prediction of the gluon-fusion cross 
section at NNLO requires the convolution with NNLO parton densities, which are not yet 
available. It is thus impossible to predict the Higgs production cross sections with NNLO 
accuracy until NNLO structure functions are accessible. 

The scale dependence at NNLO develops a similar picture as in the SM case. For large 
Higgs masses a broad maximum appears near the natural scale fi = M = M$ indicating 
an important theoretical improvement in the prediction of the Higgs production cross 
section 



3.2.2 Vector boson fusion: qq qqV*V* qqh/qqH 



h,H 



Figure 51: Diagram contributing to qq qqV*V* — > qqh/qqH at lowest order. 



Due to the absence of vector boson couplings to pseudoscalar Higgs particles A, only the 
scalar Higgs bosons h, H can be produced via the vector-boson-fusion mechanism at tree 
level [see Fig. |51[. However, these processes are suppressed with respect to the SM cross 
section due to the MSSM couplings [qv^ = sm{a — /?)/ cos(a — /?)], 

a{pp ^ qq^ qqh/qqH) = (g^^^ cr{pp ^ qq ^ qqHsu) ■ (169) 

It turns out that the vector-boson-fusion mechanism is unimportant in the MSSM, because 
for large heavy scalar Higgs masses Mh, the MSSM couplings gy are very small. The 
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relative QCD corrections are the same as for the SM Higgs particle [Fig. 27] and thus 
small fl^ . 



3.2.3 Higgs-strahlung: qq ^ V* ^ Vh/VH 




KH 



Figure 52: Diagram contributing to qq —>■ V* ^ Vh/VH at lowest order. 



For the same reasons as in the vector-boson-fusion mechanism case, the Higgs-strahlung 
off W, Z bosons, qq ^ V* Vh/VH (V = W, Z) [see Fig. is unimportant for the 
scalar MSSM Higgs particles h, H. The cross sections can be easily related to the SM 
cross sections, 

a{pp Vh/VH) = [g'i!")\{pp ^ VHsm) • (170) 

Pseudoscalar couplings to intermediate vector bosons are absent so that pseudoscalar 
Higgs particles cannot be produced at tree level in this channel. The relative QCD 
corrections are the same as in the SM case, see Fig. ^ and thus of moderate size [^ . 



3.2.4 Higgs bremsstrahlung off top and bottom quarks 




Figure 53: Typical diagrams contributing to qq/gg — > ^QQ {Q = t,b) at lowest order. 
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The scalar Higgs cross sections for Higgs bremsstrahlung off heavy quarks Q can simply 
be related to the SM case: 

a{vv ^ hQQ/HQQ) = {gi/"')' a{pp HsmQQ) (171) 



The expressions for the pseudoscalar Higgs boson ||112|| are similarly involved as the scalar 
case and will not be presented here. 

The top quark coupling to MSSM Higgs bosons is suppressed with respect to the SM 
for tg/3 > 1. Therefore Higgs bremsstrahlung off top quarks pp ^ti is less important 
for MSSM Higgs particles. On the other hand Higgs bremsstrahlung off bottom quarks 
pp — >■ Hbb will be the dominant Higgs production channel for large tg(3 due to the strongly 
enhanced bottom quark Yukawa couplings [jl8[. The QCD corrections to HQQ production 
are still unknown. 

3.2.5 Cross sections for Higgs boson production at the LHC 



Previous studies of MSSM Higgs boson production at the LHC | ]lltj|| were based on lowest- 



order cross sections or included a part of the QCD corrections. We have updated these 
analyses by including all known QCD corrections to the production processes and the 
two-loop corrections to the MSSM Higgs sector, thus rendering the results more accurate 
and reliable than in the previous studies. 

The cross sections of the various MSSM Higgs production mechanisms at the LHC 
are shown in Figs. [5^a-d for two representative values of tgjS = 1.5,30 as a function 
of the corresponding Higgs mass. The total cm. energy has been chosen as ^/s = 14 
TeV, the CTEQ4M parton densities have been adopted with as{Mz) = 0.116, and the 
top and bottom masses have been set to Mt = 175 GeV and Mb = 5 GeV. For the 
Higgs bremsstrahlung off t, b quarks, pp ^QQ + X, we have used the leading order 
CTEQ4L parton densities, because the NLO QCD corrections are unknown. Thus the 
consistent evaluation of this cross section requires LO parton densities and strong coupling. 
The latter is normalized as as{Mz) = 0.132 at lowest order. For small and moderate 
values of tg/3 < 10 the gluon-fusion cross section provides the dominant production cross 
section for the entire Higgs mass region up to M<i> ~ 1 TeV. However, for large tg/3, 
Higgs bremsstrahlung off bottom quarks, pp ^bb + X, dominates over the gluon-fusion 
mechanism through the strongly enhanced bottom Yukawa couplings. 

The MSSM Higgs search at the LHC will be more involved than the SM Higgs search. 
The basic features can be summarized as follows. 

(i) For large pseudoscalar Higgs masses Ma ^ 200 GeV the light scalar Higgs boson h 

can only be found via its photonic decay mode h —>■ 77. In a significant part of 
this MSSM parameter region, especially for moderate values of tg/3, no other MSSM 
Higgs particle can be discovered. Because of the decoupling limit for large M^ the 
MSSM cannot be distinguished from the SM in this mass range. 

(ii) For small values of tg/5 < 3 and pseudoscalar Higgs masses between about 200 and 
350 GeV, the heavy scalar Higgs boson can be searched for in the 'gold-plated' 
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Figure 54: Neutral MSSM Higgs production cross sections at the LHC = 14 TeV] 
for gluon fusion gg $, vector-boson fusion qq —* qqVV qqh/qqH, vector-boson 
bremsstrahlung qq ^ V* ^ hV/HV and the associated production gg,qq — >■ ^bb/^tt in- 
cluding all known QCD corrections, (a) h,H production for tgf3 = 1.5, (b) h,H production 
for tg(3 = 30, (c) A production for tg(3 = 1.5, (d) A production for tg(3 = 30. 
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channel H ZZ 4/^. Otherwise this 'gold-plated' signal does not play any role 
in the MSSM. However, the detectable MSSM parameter region hardly exceeds the 
anticipated exclusion limits of the LEP2 experiments. 

(iii) For large and moderate values of tg(3 > 3 the decays H,A^ t^t^ become visible 
at the LHC Thus this decay mode plays a significant role for the MSSM in contrast 
to the SM. Moreover, it will also be detectable for small values of tgjS > 1-2 and 
Ma < 200 GeV. 

(iv) For tg/3 < 4 and 150 GeV < Ma ^ 400 GeV the heavy scalar Higgs particle can 
be detected via its decay mode H ^ hh ^ 6677. However, the MSSM parameter 
range for this signature is very limited. 

(v) For tg/3 < 3-5 and 50 GeV < Ma ^ 350 GeV the pseudoscalar decay mode A — >• 
Zh l~^l~bb will be visible, but hardly exceeds the exclusion limits from LEP2. 

(vi) For pseudoscalar Higgs masses Ma ^ 100 GeV charged Higgs bosons, produced from 
top quark decays t — > H^b, can be discovered via its decay mode —>■ r^Ur- 

The final picture exhibits a difficult region for the MSSM Higgs search at the LHC. 
For tg/3 ~ 5 and Ma ~ 150 GeV the full luminosity and the full data sample of both 
the ATLAS and CMS experiments at the LHC, are needed to cover the problematic 



parameter region ||114|| , see Fig. |5^. On the other hand, if no excess of Higgs events above 
the SM background processes beyond 2 standard deviations will be found, the MSSM 
Higgs bosons can easily be excluded at 95% C.L. 



4 Summary 

In this review the decay widths and branching ratios of SM and MSSM Higgs bosons 
have been updated. All relevant higher order corrections, which are dominated by QCD 
corrections, have been taken into account. We have thus presented the branching ratios 
and decay widths of SM and MSSM Higgs particles with the best available theoretical 
accuracy. 

At the LHC the SM Higgs particle will be produced predominantly by gluon fusion 
gg —>■ H, followed by vector-boson fusion VV H {V = W, Z) and, to a lesser extent, 
Higgs-strahlung off vector bosons, V* — > VH, and top quarks, gg/qq —>■ tiH. The cross 
sections of these production channels have been updated by including all known QCD 
corrections, which are important in particular for the dominant gluon-fusion mechanism. 
Thus the final results of this review may serve as a benchmark of the theoretical predictions 
for SM Higgs boson production at the LHC. 

For Higgs masses M^ > 140 GeV the SM Higgs search at the LHC will proceed 
via the 'gold-plated' H ZZ^*^ 4/^ decay mode with small SM backgrounds. The 
extraction of four charged lepton signals at the LHC will probe Higgs masses up to about 
800 GeV. In the Higgs mass range 155 GeV < Mh ^ 180 GeV the SM Higgs boson can 
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Figure 55: MSSM parameter space including the contours of the various Higgs decay 
modes, which will he visible at the LHC after reaching the anticipated integrated luminosity 
J Ldt = 3 X 10^p6~^ and combining the experimental data of both LHC experiments, 
ATLAS and CMS [taken from Ref. M//. 
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also easily be found via its decay channel H WW I'^l'uu, by means of the specific 
angular correlations among the charged leptons. Moreover, the decay channel H l^l'vv 
may allow for an extension of the Higgs search up to Higgs masses beyond 1 TeV. For 
My/ < Mh ^ 140 GeV, the only promising decay mode seems to be provided by the 
photonic decay H 77, the detection of which, however, requires excellent energetic and 
geometric resolutions of the detectors in order to suppress the large QCD backgrounds. 
Higgs-strahlung off vector bosons or top quarks, with the Higgs decaying into a photon 
pair, may allow a further reduction of the background, but unfortunately the signal rates 
are small. In the case of excellent 6-tagging the dominant bb decay mode of the Higgs 
might be detectable, if the Higgs particle is produced in association with a W boson or 
tt pair. 

In the MSSM the neutral Higgs bosons will mainly be produced via gluon fusion 
gg — >■ $. However, through the enhanced b quark couplings, Higgs bremsstrahlung off b 
quarks, gg/qq bb^, will dominate for large tg/3. All other Higgs production mechanisms, 
i.e. vector-boson fusion and Higgs-strahlung off vector bosons or tt pairs, will be less 
important than in the SM. 

The 'gold-plated' H —>■ ZZ 41^ signal does not play an important role in the 
MSSM. On the other hand the r lepton pair decays H,A—>- t^t~ will be visible for large 
values of tg/3. The light scalar Higgs particle will only be visible via its photonic decay 
mode h 77, the branching ratio of which will be smaller than in the SM because of 
the suppressed SUSY couplings. Moreover, the decay modes H hh bb^yy, H/A 
ti and A —>■ Zh —>■ l^l^bb will be detectable in very restricted regions of the MSSM 
parameter space. Finally charged Higgs particles may be looked for in the top quark 
decays t H^b at the lower end of the charged Higgs mass range. In the search for the 
MSSM Higgs particles at the LHC, the maximal anticipated integrated luminosity will be 
needed, especially to cover the difficult region around Ma ~ 150 GeV and tg/3 ~ 5. 
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